UWTHPh-2010-4 



Emergent Geometry and Gravity from Matrix Models: 

an Introduction 

Harold Steinacker^ 

Fakultdt fur Physik, Universitdt Wien 
Boltzmanngasse 5, A-1090 Wien, Austria 



Abstract 

A introductory review to emergent noncommutative gravity within Yang- Mills 
Matrix models is presented. Space-time is described as a noncommutative 
brane solution of the matrix model, i.e. as submanifold of MP . Fields and 
matter on the brane arise as fluctuations of the bosonic resp. fermionic matri- 
ces around such a background, and couple to an effective metric interpreted in 
terms of gravity. Suitable tools are provided for the description of the effective 
geometry in the semi-classical limit. The relation to noncommutative gauge 
theory and the role of UV/IR mixing is explained. Several types of geome- 
tries are identified, in particular "harmonic" and "Einstein" type of solutions. 
The physics of the harmonic branch is discussed in some detail, emphasizing 
the non-standard role of vacuum energy. This may provide new approach to 
some of the big puzzles in this context. The IKKT model with D = 10 and 
close relatives are singled out as promising candidates for a quantum theory 
of fundamental interactions including gravity. 



harold . st einacker ©univie . ac . at 



Contents 



1 Introduction [T] 

2 Yang-Mills matrix models and quantum spaces |4] 

2.1 Quantized Poisson manifolds E] 

2.2 Moyal-Weyl solution and noncommutative gauge theory I |8] 

2.2.1 Euclidean versus Minkowski signature [TT] 

3 Geometry and gravity from matrix models [T2 

3.1 Emergent geometry [12 

3.1.1 Covariant derivatives and equations of motion [14 

3.2 Effective action, degrees of freedom and variational principle [16 

3.3 SU{n) gauge fields coupled to gravity [20 

3.4 Fermions 22 

3.5 Symmetries and conservation laws [23 

3.5.1 Tangential conservation law [2S 

4 Special geometries and perturbations 1261 

4.1 Self-duality and G^^ = g^i, [23 

4.2 Metric perturbations from U{1) gauge fields [2H] 

5 Quantization [32] 

5.1 UV/IR mixing in noncommutative gauge theory [321 

5.2 Induced gravity [34] 

5.2.1 Higher-order terms in the matrix model [36] 



6 Macroscopic description: Harmonic branch and Einstein branch 1361 

7 Harmonic branch: gravity from branes with tension 1391 

7.1 Cosmological solutions [39] 

7.2 Perturbations and Newtonian gravity from brane tension [32] 

7.2.1 Harmonic gravity bags [32] 

7.2.2 Perturbed gravity bags, Newtonian gravity and Poisson equations . . [33 

7.3 Beyond Newtonian gravity [33 

7.4 Local linearized effective action [33 

8 Conclusion [49] 



1 Introduction 

This article is intended as a pedagogical and mostly self-contained introduction to emergent 
geometry and gravity within Yang-Mills matrix models. The aim of this line of research is 
to describe the physics and the geometry of general noncommutative (NC) spaces in these 
matrix-models, and to assess their viability as a quantum theory of space-time including 
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matter, gauge fields and gravity. We hope that the present review provides a useful basis for 
further research in this direction. 

Our starting point is the identification of a gravity sector within noncommutative gauge 
theory. NC gauge theory has been considered previously as a deformation of Yang-Mills 
gauge theory, living on NC space. From that point of view, it is well-known that the U{1) 
sector of U{n) gauge theory on the Moyal-Weyl quantum plane (which is the simplest 
example of a NC space) plays a special role: it does not decouple from the remaining SU{n) 
degrees of freedom, and its quantum effective action is drastically different from its commuta- 
tive counterpart due to UV/IR mixing. These and other "strange" features have been viewed 
as obstacles for the physical application of NC gauge theory, and a relation to gravity has 
been widely conjectured. This conjecture is corroborated further through string-theoretical 
matrix models such as the IKKT model [l], the description of gravitons in the BFSS matrix 
model and related models (2|-|6|, and matrix-models describing curved or "fuzzy" quantum 
spaces (3]10)[12)[13|. 



More recently, Rivelles observed 14 that the coupling of matter to NC U{1) gauge fields 



can be rewritten in terms of an effective metric, thus pointing out a direct relation between 



NC gauge theory with geometry and gravity. Such a relation was also advocated in 15, 16 



in particular a relation between self-dual gauge fields and gravitational instantons. This 



motivated a series of papers 17 22 studying the effective geometry of general NC branes 



of Yang-Mills matrix models, which provide a non-perturbative formulation of NC gauge 
theories. The central point in these works is that NC gauge theory is in fact a theory of 
noncommutative spaces and their fluctuations. The geometry and the effective metric is 
not a fundamental degree of freedom which is put in by hand, rather it arises effectively in 
the low-energy or semi-classical description of the model. That is the reason for calling it 
"emergent NC gravity". 

In order to have a well-defined framework, we will focus on matrix models of Yang-Mills 
type. These models have non-commutative spaces or space-time as solutions, i.e. quantized 
Poisson manifold. Thus space-time and geometry are dynamical rather than put in by 
hand, and the models should be considered as background independent. U{1) fluctuations 
of the matrices around NC space-time describe geometrical deformations such as gravitons, 
while SU (n) fluctuations describe nonabelian gauge fields. The kinetic terms of these fields 
arises from the commutators in the matrix model, and encodes an effective metric which 
is essentially universal for all fields and matter. Since this metric is dynamical, it must be 
interpreted in terms of gravity. This leads to an intrinsically non-commutative mechanism 
for gravity, combining the metric and the Poisson structure in a specific way. It provides a 
natural role for non-commutative or quantized space-time in physics. 

Let us describe some results in more detail. Space-time is described as a 3+1-dimensional 
NC brane Aie C M.^ (possibly with compactified extra dimensions), which carries a Poisson 
tensor 6^'^{x). All matter and gauge fields live on this space-time brane, and there are no 
physical fields propagating in the ambient D-dimensional space unlike in other braneworld 
scenarios such as [23]. An effective metric ~ ^^'^ 9^" Qij^iyi arises on this space-time brane, 
which governs the kinetic term of all fields more-or-less as in general relativity (GR). This 
metric is dynamical, however it is not a fundamental degree of freedom: it is determined 
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by the embedding A4g C M'^, and the Poisson tensor 6'^'^ describing noncommutativity. 
Hence the fundamental degrees of freedom are different from GR, and can be interpreted 
alternatively in terms of NC gauge theory. This makes the dynamics of emergent NC gravity 
somewhat difficult to disentangle, and the effective metric is not governed in general by the 
Einstein equations. 

We will identify 2 classes of solutions: in the "Einstein branch", solutions of the Einstein 
equations can be realized as embedded submanifolds for D > 10. Since the Einstein-Hilbert 
action arises upon quantization and is not part of the bare matrix model, the model must 
be free of UV/IR mixing above a scale A identified as Planck scale. This singles out the 
IKKT model [I] or close relatives, with D = 10 and maximal supersymmetry above A. In 
contrast, the solutions in the "harmonic branch" are governed by the brane tension rather 
than the induced Einstein-Hilbert term. That branch is interesting because the physics of 
vacuum energy is different from GR, providing a new perspective notably for the cosmological 
constant. For example, a cosmological solution is discussed which requires far less fine-tuning 
than the ACDM concordance model, featuring a big bounce. Newtonian gravity arises in 
this branch through an interesting mechanism, with long-distance modifications which may 
be relevant for the issue of dark matter. Even though the solar system precision tests appear 
to be violated in the present ansatz, it seems conceivable that some of these features carry 
over into a modified Einstein branch. 

Apart from possible deviations from GR and its intrinsic appeal, the main motivation for 
the matrix model approach is that notably the IKKT model promises to provide an accessible 
quantum theory of gravity as well as the other fundamental interactions. The reason is that 
it can be equivalently interpreted as TV = 4 NC SYM on Mg. Thus it is (expected to be) UV 



finite, and free of pathological UV/IR mixing 24,25 . The remaining mild UV/IR mixing 
below the TV = 4 SUSY breaking scale is nothing but induced gravity and therefore welcome. 
In particular, the quantization amounts to the quantization of non-commutative Yang-Mills 
action, rather than the Einstein-Hilbert action. This elegantly circumvents many of the 
technical problems in quantizing GR. Hence these models under consideration define some 
accessible quantum theory of gravity, which is at least not too far from real gravity. 

Studying the geometrical properties of the matrix models, we will obtain structures and 
features which are familiar from string theory. For example, the effective metric on the brane 
is very similar to the open-string metric on branes in analogous string theory considerations. 
This is to be expected in view of the relation with string theory. However, the matrix model 
is taken as a starting point here, because it is intrinsically non-perturbative and extremely 
simple. Hence we cannot just use results from string theory, rather everything will be derived 
in a self-contained manner starting from the matrix model. 

Some words on the "philosophy" of emergent NC gravity are in order her^ Rather 
than trying to carry over some formulation of GR into the context of quantized space^ one 
starts with a model whose fundamental degrees of freedom are different from GR, but which 
provides a physically viable description of geometry and gravity in some geometrical phase 



26 



^For a review of analogue models for emergent gravity see 
^which would not simplify the problem of quantizing gravity in any obvious way. 
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or limit. Matrix models are not only very natural in the context of quantum spaces, they 
appear to realize the idea of emergent gravity, and moreover - most importantly - promise 
to define a good quantum theory. 

In order to provide a coherent picture, we focus on the matrix model approach following 



17-20 and provide the necessary background, rather than attempting a survey of related 



work and other points of view. In particular, the relation with string theory is discussed only 
briefly in section [3| This allows to keep the presentation simple and accessible, but does not 
imply any disregard of other results which are not covered. The reader is assumed to have a 
basic background in high-energy theory, gauge theory and general relativity, but no previous 
knowledge in noncommutative field theory is required. 

The structure of this review is as follows: in sections 2 and 3 the basic geometrical 
properties of noncommutative branes are derived. In section 4 a special class of geometries 
is discussed which appears to be most relevant for physical applications. The quantization 
(focusing on one-loop) is discussed in section 5. The equations of motion for gravity coupled 
to matter are then discussed in section 6. Section 7 is more speculative in nature, exploring 
the physics of the "harmonic branch" of solutions which is somewhat different from GR; 
section 7.4 contains some new material. The exploration of the more conventional "Einstein 
branch" and possible other scenarios is left for future work. 



2 Yang-Mills matrix models and quantum spaces 

We consider the following type of Yang-Mills matrix models 

S = -^^Tt{[X^,X'][X^' ,X'']r^aa'Vw) (1) 

where 

r/a„/ = diag(-l, !,...,!); (2) 
the Euclidean version of the model is obtained by replacing 

Vaa' Saa'- (3) 

The degrees of freedom of this model are hermitiar]|^ matrices X° G Mat(oo,C) for a = 
0, 1,2, ...,D — 1. We introduced also an energy scale Aq which gives the matrices X" the 
dimension of length. The action is invariant under the fundamental gauge symmetry 

X^" UX^'U-\ U e W(oo) (4) 

as well as a global rotational and translational symmetry, where rotations act as X" — > X* 
for A^ G SO{D - 1, 1) resp. SO{D), and translations X" X" + c"! for G M. However 
there is no space-time or geometry to start with; space and geometry only emerge on a 

■^in the Minkowski case, we will assume that the time-like matrices are anti-hermitian. This will be 
addressed in more detail below. 
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given solution of the model, and are therefore dynamical. The only geometrical input is 
the constant metric rjab resp. 6ab for the "embedding space" R^. This space is however 
unphysical and only serves to "suspend" the physical space-time brane and its emergent 



geometry, as we will see. Fermionic matter will be added in section 3.4 

This class of models can be obtained as dimensional reduction of Yang-Mills gauge theory 
to a point. In particular, the IKKT model with D = 10 is singled out by an extended ma- 
trix supersymmetry [l]. This is probably essential to obtain a well-defined quantum theory. 
Some modifications of these models are conceivable, such as a "mass" term Tr{m?riabX°'X^), 
or cubic terms Tr{Q ahcX^'iX^ ■, X'^]) ■ These "soft" terms will have mild impact on the quan- 
tization properties on a 4-dimensional space-time brane, and allow to obtain interesting 
compactification scenarios. 

The equations of motion corresponding to the action ([T]) is obtained by varying the entries 
of the matrices X", which gives 

[X^[X^X^']]r/,,, = 0. (5) 

This equation has various types of solutions. An obvious solution is given by any commuting 
set of matrices = 0. These "commutative" configurations do not support propagat- 

ing fields in this model and are in some sense singular; we will not pursue them any further 
here. The prototype of the class of solutions which is important here is given bjj^ 

X«^f^n = f^M, /i = 0,...,2n-l, 2 = l,...,D-2n (6) 



where are generators of the Moyal-Weyl quantum plane Mg", 

[X^',X''] = ie^"'X fi,u = 0,...,2n-l. (7) 

Here 6^'^ is a constant antisymmetric numerical tensor. This solution should be interpreted as 
quantization of fiat M^" C M.^ with Poisson structure {x'^,x'^} = 6^'^. From a mathematical 
point of view this is very similar to quantum mechanics, where phase space is quantized using 
the same commutation relations. The matrices X^ should thus be interpreted as quantized 
coordinate functions on M^". The matrices 0* can similarly be interpreted as quantized 
functions resp. scalar fields on M^"^, which happen to vanish in this particular solution. 
As in quantum mechanics, the above commutation relations imply space-time coordinate 
uncertainty relations 

This means that space-time becomes "fuzzy" at the noncommutative scale defined by 

A^"^ = (detr'')-^/^ (9) 
Such uncertainty relations can be motivated by standard arguments combining quantum 



gravity with mechanics 27 



The important message is that space-time arises as a solution of the matrix model, it is 
not a fixed background. This suggests some relation with gravity. 

''^The particular embedding of M^" C implied by this choice of variables is of course arbitrary, and 
will be generalized below. 
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2.1 Quantized Poisson manifolds. 

To proceed, we recall the concept of the quantization of a Poisson manifold (A^, {.,.}), 
referring e.g. to ^28j and references therein for more mathematical background. A Poisson 
structure is an anti-symmetric bracket {., .} : C{Ai) x C(A^) — )■ C{A4) which is a derivation 
in each argument and satisfies the Jacobi identity, 

{fg, h} = f{g, h} + g{f, h}, {/, {g, h}} + cycl. = 0. (10) 

It is sometimes useful to introduce an expansion parameter of dimension length"^ and write 

{x>',x''} = 9'"' = ee^''ix) (11) 

where Oq'^{x) is some fixed Poisson structure. Given a Poisson manifold, we denote as 
quantization map an isomorphism of vector spaces 

X: C{M) ^ Me C. Mat(oo,C) . . 

fix) ^ F ^^^^ 

which depends on the Poisson structure I ^Xg, and satisfied 

Xifg)-Xif)Xig) and ^ (x(z{/, (7}) - [X(/),X(^?)]) ^ as ^ 0. (13) 

Here C{Ai) denotes some space of functions on A^, and A^g is interpreted as quantized 
algebra of function^ on Ai. Such a quantization map X is not unique, i.e. the higher-order 
terms in (13) are not unique. For example, the Moyal-Weyl quantum plane M.f^ ([t]) is a 
quantization of M^" with constant Poisson tensor 6'^'^, and a natural ("Weyl") quantization 
map is given by 

X: C(M2«) ^ M2n ^ Mat(cx),C) 

The map X allows to define a "star" product on C{Ai) as the pull-back of the algebra resp. 
matrix product in Aig, 

f ^ g := X-\X{f)X{g)) (15) 

which allows to work with classical functions, hiding 9 in the star product. In the example of 
Mg", it leads to the well-known Groenewald-Moyal star product on M^". The compatibility 
(13) with the Poisson structure is encoded by the requirement 

f^9 = f9+\{f,g} + 0{6^). (16) 

Kontsevich has shown [29] that such a quantization always exists in the sense of formal 
power series in 9. This is a bit too weak for the present context since we deal with operator 



^The precise definition of this limiting process is non-trivial, and there are various definitions and ap- 
proaches. Here we simply assume that the limit and the expansion in 9 exist in some appropriate sense. 

is the algebra generated by = I(x^), or some subalgebra corresponding to well-behaved functions. 
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resp. matrix quantizations. However we will assume that 6^'^ is non-degenerate defining a 
symplectic structure 



1 
2 

for which stronger mathematical results are available. 



oj = ^e-^dx'' A dx", (17) 



De-quantization and semi-classical limit. Since the matrix model provides directly 
quantum spaces such as Mg", we are actually faced with the opposite problem of "de- 
quantization", i.e. the semi-classical limit of a quantum space. The above concepts can be 
used as a tool to understand solutions of the matrix model in the language of ordinary spaces 
and effective actions. Denoting the inverse of the quantization map as fg : = Xg^{F), we can 
replace [F,G] by [fe,g9]-k with X being understood. Now we can define the "semi-classical 
limit" as the leading non-vanishing term in an expansion in 6, dropping all higher-order 
terms in 6; this semi-classical limit will be indicated by ~. In particular, 

[F,G] ~ t{f,g} i+0{9')) (18) 

which will be used throughout. Thus one can simply replace commutators with Poisson 
brackets in the semi-classical limit, and any ambiguity in Xq would show up only in higher- 
order corrections; this is familiar from quantum mechanics. In particular, if some generators 
X^^ generate the entire operator resp. matrix algebra, we can at least formally write every 
"smooth function" on A^e as F = F{X^). Then 

~ x^ (19) 

can be interpreted as quantized coordinate function on A^, and we can write 

[X^F]~z^^^(x)^/(x). (20) 

The important point is that matrices F G Mat(oo, C) should be interpreted as quantizations 
of functions on A^, and their non-commutative product encodes an energy scale A^vc- It 
is also helpful to recall the concept of coherent states or optimally localized states, which 
saturate the uncertainty relations ([s]). They make explicit the fact that on length scales larger 
than L^c = ^~nci quantized spaces Aig look and behave like their (semi-) classical limit 
M.. In much of this paper we will restrict ourselves to the semi-classical limit of the matrix 
model and its geometrical solutions. 

Finally, under favorable assumptions the integral on Al is related to the trace as follows 
(cf. |28|) 

= /»(iF^<^' ~ '''' 

e{x) := Pfaff(^^"(x))^ = |0^-(x)|2^i =: L]^c{x). (22) 
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The symplectic volume form is singled out by the compatibility condition 

Tr[X(/),X(^)]~ J ^{L9} = 0. (23) 

It is instructive to recall here the Darboux theorem, which states that every symplectic 2- 
form is locally constant in suitable coordinates. This means that in suitable coordinates, any 
quantum space (with non-degenerate 6'^'^) locally looks like Mg". It also provides an intuitive 
justification for ([2T|, which then reduces to the Bohr-Sommerfeld quantization rule. 



2.2 Moyal-Weyl solution and noncommutative gauge theory I 

Let us start with the Moyal-Weyl solution ([t]), and recall that translations on M^" are realized 
as inner automorphisms. 



Uf{X^)U-' = fix" - e^'^k^), U = e'"^"" . (24) 
This justifies to define a partial derivative operator on Mg" in terms of a commutator 

d 



5;^,/ = -^^."[^^/]' (25) 

consistentj^ with the semi-classical limit (20). The fact that derivatives are inner deriva- 
tions constitutes the basic difference between the commutative and the non-commutative 
framework. This is compatible with the definition of the integral over as 

/ /^=^Tr/, = ^\e^'^{x)\ = A^t; (26) 

which has the correct classical limit (21), and satisfies Jj^4 d^f = up to "boundary terms". 

Now consider fiuctuations of the above solution M.j C ([6|, ([t]) of the matrix model, 
focusing on the 4-dimensional case to be specific. Such fiuctuations can be parametrized as 







(27) 



where A'^ and 0* are "small" hermitian matrices. As explained above, they can be interpreted 
as smooth functions on M^, i.e. A^ = A^{X) ~ A^{x) and 0' = 0*(X) ~ 0*(x). Consider 
now the change of variables 

Af" = -e'^'A^ (28) 



where A^, is hermitian. Using (25) one finds 



= [X^ + A^,f]=zr''{^J + t[A,,f]) ^ te^^DJ, 



— ,,,,/ — , ,,,/ , — 



[x^xn = -^r^r''(r,t, + F^v) (29) 



''the Poisson structure is always assumed to be non-degenerate. 
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where F^y = d^A^ — d^Ai, + i[Af^,A,^] is the U{1) field strength on Mg. The symmetry ^ 
of the matrix model now acts on the fiuctuation fields as 

A^, UA^U~^ + Udf,U-\ 0* ^ U(f)'U-\ (30) 

This clearly has the structure of gauge transformations of Yang-Mills gauge field^ Moreover, 
the action ([T]) can now be rewritten as 

S = S[A] + = / d^x ^ (&^' G'"'' F^v + G'^'^v^u 

+lG'^''D,<pWy^^,j + I (31) 
dropping surface terms, where 

Qi^y = J-r'^'^'^'r^ IG'^^I = 1. (32) 

These formulas are exact (up to boundary terms) and define some noncommutative field 
theory 31 . Thus the matrix model ([T]) on the solution can be viewed as a noncommu- 



tative gauge theory coupled to scalar fields (jf in the adjoint, with effective flat metric G^^ . 
Redeflning (p^ = 0* gives the scalar flelds the usual dimension of mass. 

Nonabelian gauge fields. Nonabelian gauge flelds arise very naturally in the matrix 
model. Observe flrst that for every solution X" of 

X'^ ® 1„ (33) 

is also a solution, naturally interpreted as n coincident branes. The fluctuations around such 
a background can be parametrized as 

(34) 

considering to be speciflc. In analogy to the above considerations, it is easy to see that 

= _r-A,,„(X)A", 0* = 0l(X)A'^ (35) 

describe u(n)- valued gauge resp. scalar flelds on M^, denoting with A" a basis of u(n). The 
matrix model action for the fluctuations then takes the form 

S = / ciV^trfG'^^'G'^'^'F^.F^v' +G'^"r?M-ll„ 

+ 1 D,<I)'D,(I)^,, + 1 [0\ 0^] [0^ 0^']5,,,5,,v) (36) 



^The can be interpreted as "covariant coordinates" ISOl. 



where tr() denotes the trace over the u(n) matrices, and F^^ is the u(n) field strength. 

This mechanism is very remarkable: It shows that gauge theory arises "automatically" 
in the matrix model, there is no need to define sophisticated mathematical structures such 
as principal fiber bundles and connections. This is one of the reasons why these models are 
so interesting. On the other hand, even though the derivation and interpretation in terms of 
gauge theory on seems impeccable, it is nevertheless misleading and physically "wrong" 
for the U{1) sector, for a number or reasons: 

• While U{n) gauge fields arise naturally on M^, there seems to be no way to separate 
them into U{1) and SU{n) gauge fields. This is fundamentally tied to the noncommu- 
tativity of space. For example, if we were to impose the constraint that is traceless 
tiA^ = 0, then a gauge transformation 6A^ = (9^A + [A^, A] will re-introduce a trace- 
U{1) component ~ !„. Indeed the trace-t/(l) components of gauge fields and scalar 
fields are in-separably entangled with all other fields. We will see later that this is 
nothing but a gravitational coupling, and the U{1) components will be understood as 
geometrical degrees of freedom. Only the SU{n) components of the gauge field turn 
out to be physical gauge fields, while the trace-?7(l) components turn into gravitational 
waves. The latter will be denoted as "would-be U{1) gauge fields" henceforth. 

• We will see in section |5] that upon quantization, the one-loop effective action contains 
"strange" UV/IR mixing terms for the trace-[/(l) component, but not for the SU{n) 
terms. This indicates that the would-be U{1) gauge fields should not be interpreted as 
photons. These "strange" terms will be understood as induced gravitational terms. A 



notable exception is the = 4 SUSY case, which contains no such terms 25 



• Eq. (24) says that translations on Mg are nothing but particular gauge transformations. 
Thus the U{1) gauge transformations on noncommutative space contain space-time 
transformations, which can be interpreted as symplectomorphisms which form a sub- 
group of the volume-preserving diffeomorphisms. In the same vein, there are no local 
observables in noncommutative gauge theory: for example, the U{1) field strength F^^^, 
is not invariant under U{1), but transforms as a scalar field under these symplectomor- 
phisms. Gauge invariant observables such as Tr(e**^'""^'') involve an integral over space. 
All these are signatures of gravity. 

• Slightly modified versions of the matrix model have quantum spaces with non-trivial 
geometries such as fuzzy spheres as solutions, cf. [7||To|[lT|[32]. It is then obvious that 
fiuctuations of these solutions, i.e. fiuctuations of the matrices, should correspond to 
fiuctuations of the geometry rather than U{1) gauge resp. scalar fields. 

The resolution of this puzzle is the key for identifying emergent gravity in the matrix 
model. In the geometrical interpretation explained below, the trace- U{1) components of 
Af^{x) and will completely absorbed in the effective metric G^'^{x), leaving only physical 
SU {n) gauge fields and scalars as well as fermionic matter coupled to gravity. 
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2.2.1 Euclidean versus Minkowski signature. 



Let us take a closer look at the effective metric G^*^ for gauge fields on M^, focusing first on 
the Euclidean case. Using a SO (4) rotation, we can assume that 6'^'^ has the canonical form 



f ° 








—a 


\ 




























y a 











/ 



(37) 



Clearly the corresponding symplectic 2-form u = 9~^dx^ A dx'^ is (anti-) self-dual -ku = ±u 
if and only if = 1, where * denotes the Hodge star defined by e'^'^p"' and on M^. On 
the other hand, (32) gives 

G-M- = 1 ^''^'r^' V.' = diag(a^ a 



-2 -2 2\ 



(38) 



Therefore G^'^ = 5'^^ if and only if u is (anti-)self-dual, otherwise G'^'^ differs from the 
embedding metric 6^'^ via a volume-preserving linear transformation. 

Now consider the physical case of Minkowski signature with embedding metric r/^^ = 
(— 1, 1, 1, 1). The standard point of view in NC field theory is that 6'^'^ is real-valued. Using 
suitable Lorentz transformations, we can again assume that it has the above form (37), so 
that 

G^"" = 4 O^^^'d^'^'r]^,^, = dmg{a\ a'^, a'^, -a^). (39) 



9^ 



This indeed has Minkowski signature, but the role of time has switched from the first to the 
last coordinate. In particular there is no way that the effective metric G^'^ agrees with the 
embedding metric r/^^. This is a priori not a problem since only G'^'^ couples to the physical 
fields. However it is quite counterintuitive, and becomes more problematic at the one-loop 
level as we will see. To gain some insight consider again the Euclidean case and write 



with real . Having in mind a Wick rotation of the matrices 



(40) 



(41) 



such that the coordinate functions are related as usual via x° — )■ it, the above commutator 
becomes 

[T, X^\ = 6% =: te%, [X\ X^] = id% =: tff^. (42) 

This leads to a generalized Poisson tensor with imaginary time-like components, accord- 
ing to (X'^)* = r]^yX" . We can again assume that it has the form 



/ 











—ia 


\ 
































\ 


ia 











/ 



(43) 
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in suitable coordinates. Now the effective metric (32) has the form 



G'"" = I2 «V^' = diag(-«2, a-', a'', a'), (44) 

and the time coordinate is in the expected slot. Moreover, -kcu = ±a; if and only if = 1, 
where -k denotes the Hodge star defined by rj^^^ on M.^; here we adopt the convention that the 
Wick-rotated epsilon-tensor is given by e'^^^s _ ^_ Thus we again have G'^'^ = rj^^ if and only 
if u is (anti-) self-dual -ku = iw, as in the Euclidean case. This seems to be the appropriate 
concept of Wick rotation in the matrix model framework, which will generalize to the case 
of non-trivial geometries and gravity. 



3 Geometry and gravity from matrix models 

We have seen that NC space-time is obtained as a solution of the matrix model. Hence space- 
time is dynamical, which strongly points at gravity. On the other hand, the fluctuations of 
space-time were identified as gauge fields on M^. The key observation of emergent NC 
gravity is that these apparently different interpretations are two ways of looking at the same 
thing. This leads to an intrinsically non-commutative and compelling mechanism for gravity. 



realizing ideas in 14 ,15 . Here we follow a systematic approach and derive everything from 
the framework of matrix models. 

We consider configurations X" of the matrix model which can be interpreted as quantized 
embedding functions 

X'^ ~ : Mf" ^ (45) 
determining some 2n-dimensional submanifold ("brane") of M^. Here ~ denotes the semi- 



classical limit as discussed in section 2.1 The intuitive picture is that there are sufficiently 
localized ("coherent") states such that (X") spans the manifold A^^" C M^. We can then 
split (at least "locally") the set of matrices as in ([6]), 

X" = ( ] , /i = 0,...,2n- 1, i = l,...,D-2n (46) 



such that the 2n generators generate the full matrix algebra[j Mat (00, C), and therefore 
0* = (j)^{X^) ~ (j)^{x^). Since the commutators always satisfy the Jacobi identity, 

they can be interpreted as quantizations of a Poisson structure on Ai, 

[X", X'] =: lO^'^iX^') ~ lO^'^ix), fx,u = l, 2n (47) 

Thus Mat(oo, C) = Cei^Ai) is interpreted as quantized algebra of functions on a manifold M. 
with Poisson structure 9^^^{x). In particular, the X^ ~ x^ define quantized (local) coordinate 
functions on M., which we will denote as "matrix coordinates". 



The difference to the discussion in section 2.2 is that the "would-be U{1) gauge fields 

1-1 



are absorbed in X^, avoiding the unphysical splitting of the Poisson tensor 9 J; in (29). 



Since Mat{oo,C) is irreducible, this implies that 6^" as defined in (47) is non-degenerate. 
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Note that all physical fields in this framework arise from fiuctuations of the matrices 
around such a background (leading to nonabelian gauge fields and scalars), and from the 
fermionic matrices \E'. Since Mat{oo,C) = Cg{Ai), it follows that they all live only on 
the brane Ai, hence there is no physical higher-dimensional "bulk" which could carry any 
propagating degrees of freedomp^ Moreover all fields propagate according to an effective 
metric, which we now identify. 



3.1 Emergent geometry. 

To understand the effective geometry of A^^", consider a test-particle on A^^", modeled by 
a scalar field ip for simplicity (this could be e.g. an su{k) component of 0*). In order to 
preserve gauge invariance, the kinetic term must have the form 



S[^] 



-Tr[X^^][X^^]r/„, 



(2vr) 



1 



{2'kY 
1 

{2'kY 



denoting the D natural vector fields on Ai defined by the matrix model as 



where 18 



G^"(x) 

^-(n-l)a 



'v'[X) 



e-''d^^\x)d'''\x)g^.,^ 



9 



\9^^u[x) 



9' 



(48) 

(49) 

(50) 
(51) 

(52) 



Here g^v{x) is the metric induced on C MP via pull-back of rjab- The normalization factor 
e~" is determined uniquely such that 



1 

9^ 



G 



flU I 



(53) 



except for n = 1 which we exclude for simplicity. Therefore the kinetic term for ip on 
Ai is governed by the metric G^^{x), which depends on the Poisson tensor 9^'^{x) and the 
embedding metric g^^{x). We will see below that the same metric also governs nonabelian 
gauge fields and fermions in the matrix model (up to possible conformal factors), so that 
G^u must be interpreted as gravitational metric. There is no need and no room for invoking 



10 



This does not exclude the existence of compactified physical extra dimensions in the matrix model. 



13 



any "principles". Here dim^ = dim^'^'^ = L^, hence e ^ has the dimension L ^" set by the 
noncommutativity scale. We finally note that 



\G^y{x)\ = \g^,u{x)\, 



2n=4 



(54) 



which means that in the 4-dimensional case, the Poisson tensor 6^'^ does not enter the 
Riemannian volume at all. This is important for stabilizing fiat space as we will see, and is 
one of several reasons why 4 dimensions are special in this framework. 

3.1.1 Covariant derivatives and equations of motion. 

To understand the meaning of the matrix equations ([s]), we note that the double commutator 
[X"', [X'^, (p]]riab in the semi-classical limit reduces to the covariant Laplacian corresponding 
to the effective metric G^^y. This is based on the following facts: 

Lemma 1 Using the above definitions for the metrics g^uiG^'^ and the (non- degenerate) 
Poisson structure 6^'^ on M. c MP with Cartesian embedding coordinates : — M, a = 
1,...,D, the following identities hold: 



(55) 
(56) 



Here ip is a scalar field on M. c M''^, Vg resp. Dg denotes the covariant derivative resp. 
Laplacian corresponding to G^y, and 



yG^'^{x)g^,{x). 



Proof We first note the following useful identities for the metric G^'^: 

1 







-e"e-"9^(e"G"^^^) = -e-"r'^'9^(rX^,(a;)). 



(57) 

(58) 

(59) 



The first is a consequence of the Jacobi identity; for the short proof see Appendix A. (59) 



follows from (50) and (58). Now we compute 

{x^{x^y.}K, 



e^^pd,x''d,{e-^dyd^v)r,ab 
e^pd,{d^x''e'''^d,x'dr,^)r]ab 
e^^dpie^^g^yd,^) 

e^^e-^g^ydpd^^ + e>^pdp{e"^g^,)d^^ 



(60) 



using (58) and (59). Finally, (56) is shown in Appendix B ( |227 ), based on a result in (2l] . 



14 



Using these identities, the matrix equations of motion [X^, [X'^, X^']]rjbb' 
classical limit become 



0, 



in the semi- 

(61) 
(62) 



independent of the splitting X" = (X^,0*) into coordinates and scalar fields. Note that 



is viewed as scalar field on Ai here. Together with (62), (56) can be written as 



V^(e^ 



(63) 



These constitute the covariant equations of motion of the bare matrix model without matter. 
It is remarkable that (63) is completely "intrinsic", i.e. the embedding of M. plays no role. It 



provides the relation between the noncommutativity Q^^{x) and the metric G^^ . Since (63) 



has essentially the form of covariant Maxwell equations, it should have a unique solution for 
a given "boundary condition" such as 



(^^iu{x) — )■ 9^1, = const for |a;| — )■ oo 



(64) 



up to radiational contributions, which will be identified as gravitational waves. We will see 



that (63) is in fact a consequence of a conservation law resp. Dyson- Schwinger equation. 



and is therefore protected from quantum corrections [21]. 

In general, any 2n of the matrices X° ~ x"" can be considered as coordinate functions 
on A^; hence the matrices define preferred "matrix coordinates". (62) implies that these 

matrix coordinates satisfy the harmonic gauge condition 



0. 



(65) 



From the point of view of GR, (65) would be interpreted as gauge fixing condition, disposing 
of diffeomorphism invariance which does not make sense in the matrix model. Since gauge- 
dependent objects are always unphysical, this has no implications on the physical content of 
the model. The equations of motion in these matrix coordinates are given in 18 . 

Finally a comment on Lorentz invariance is in order. Even though 6^'^ explicitly breaks 
Lorentz invariance (just like any other background field), the physical fields of the model do 
not couple directly to 9^'^, and the effective actions such as (48), (83) are (locally) Lorentz- 
invariant with respect to G^u. There is no field in the model which is charged under the 
trace- U{1) gauge field whose flux corresponds to Therefore Lorentz-breaking effects 



can arise only through higher-order corrections in 9^'^ or through quantum effects, both of 



which are expected to be suppressed by powers of A^^, or probably A^y"^. 



^^In fact the ncxt-to-leading order corrections are expected to be suppressed by 0{6^), because one can 
choose a star product 29 such that [f,g] = i{f,g} + 0{9^). 
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3.2 Effective action, degrees of freedom and variational principle 

We now show how to obtain the above geometrical equations of motion directly from the 
semi-classical limit of the matrix model action ([T|, which using (53), (57) can be written as 



Before applying the variational principle, we should first clarify the independent degrees of 
freedom. 



Geometrical degrees of freedom. The semi-classical action (66) depends on the Poisson 
tensor 9^'^{x) and the embedding metric g^p = 7]^^ + dfj_(f)'^du(f)^6ij, which in turn depends on 
the embedding functions (f)^{x^). On the other hand, the fundamental degrees of freedom 
are the D matrices resp. functions X"", which can be varied independently. Semi-classically, 
these variations dX"" can be decomposed into tangential and normal fluctuations w.r.t. the 



background brane M.. This is most transparent in normal embedding coordinates (99), 
where 50* G TAi^ resp. SX^^ G TA4 correspond to transversal resp. tangential variations of 
the brane. Clearly the variations of the scalar fields change the embedding A4 C M.^ 

and the embedding metric g^^, while tangential variations 

{x''y = x^' + A^ (67) 
lead to a variation of the Poisson tensor on Ai according to 

e^"" ~ -i[x^, x""] -i[{x^y, {x")'] ~ e^"" + e^'^dpA'' - o'^^d^A^ - i[A^, A''] 
= {1 + A ■ d)e^"' - e^'^^'e"'"' F^,^. 

= {e"")' (68) 

using the Jacobi identity, where F^/^/ = df^'A,^/ — dyA^i and A-d = A^dp. For the symplectic 
2-form, this amounts to 

{Opiy = e-J - e-l^e-lM^'' = {i + A- d)e-J + f,.. (69) 

Note that these horizontal variation determine a vector field A = A^d^ and thus an (in- 
finitesimal) diffeomorphism — A^. In particular, the f/(oo) gauge symmetry (|4]) resp. 
its infinitesimal version X" — (X")' = X" — 2[A,X"] and 0* — )■ 0* — z[A, 0*] corresponds to 

5X^ = ~ -e^'^^A{x), 6(f)' ~ -e^"'dp(f)'d^X{x). (70) 

This defines an exact symmetry of the matrix model corresponding to (infinitesimal) sym- 
plectomorphisms Sympe 

{A,.} = (9,A)r''9^. (71) 
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By definition, symplectomorphisms leave 9^^{x) invariant C_a^uj = 0, and form a subgroup 
of tlie group of volume-preserving diffeomorpliisms due to Liouvilles theorem. From the 
point of view of U{1) gauge theory on Mg, they correspond to gauge transformations — )■ 

+ d^\{x), thus providing a connection between gauge theory and (emergent) gravity. 

A more geometric understanding is obtained as follows: since locally the generate the 
full matrix algebra, a general variation of the matrices X"'' = X"- + A"" ~ x°'{x^) + A°'{x^) 
corresponds semi-classically to a modified embedding 

X"' : M' ^ M ^ (72) 

in terms of some map ^ , i.e. (x")' = o v]/. In general \l/ will be neither an isometry nor a 
symplectomorphism, and the new configuration X"-' is completely inequivalent to X"". If \E' 
happens to be a diffeomorphism — )■ (in particular it preserves the embedding), then 
the induced metric g^y is equivalent, but the Poisson structure and therefore G^" is not. 



This is the case of tangential transformations (67). In particular, gauge transformations 
^ define such a symplectomorphism ^! on A^, hence {x" o^i^x^ o = {x°',x^} o v]/^ and 
the effective geometry is the same. The matrix model action and its semi-classical limit are 
invariant under this symmetry, because the symplectic volume resp. the trace is preserved. 
Hence in the matrix model, transformations should always be considered as "active" as 



in ( 72 ) , rather than as diffeos in some given coordinate system. All behave as scalar fields 
on M. and are on equal footing. This is why gauge transformations (|4]) (and more generally 
the transformations (67)) do not appear to act as Lie derivatives; for example, (69) is not 



the Lie derivative of a 2-form along the vector field A which would be jC^o; = F^ydx'^ A dx'^. 

The bottom line is that the 4 tangential degrees of freedom A^ are transmuted into 
the 3 (resp. 2 "on-shell") degrees of freedom of a; = \6~ldx^ A dx^ resp. the U{1) field 
strength F^^, plus one gauge degree of freedom corresponding to symplectomorphisms. This 
means that one can consider the embedding cjf and 9^^^ as independent geometrical degrees 
of freedom, and use a semi-classical variational principle for the actions of the type 

S = J d^''xpC{(l)\x),9'"'{x)) (73) 

as illustrated below. Even though this gives the correct equations of motion, one should 
keep in mind that the fundamental degrees of freedom are the matrices X"-. Their variation 
does not simply separate into 69^" and 50' but involves also a change of coordinates. 

Semi-classical derivation of the equations of motion. Once the classical degrees of 



freedom are understood, one can work with the effective action (66) for the matrix model 
in terms of a sub-manifold Ai C M"^, and derive covariant equations of motion for the 
embedding and the Poisson structure. As we just explained, the most general variations of 
the latter can be parametrized in terms of a U{1) gauge field as 

= V^M, - VJA^, (74) 

^^It is in general much smaller that the volume-preserving diffeos except in 2 dimensions, hence a com- 
parison with unimodular gravity is quite inappropriate. 
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and the variation of the embedding can be captured by 50*. Thus the variation of the 
effective action (66) is 



SS 



I J d'^x^f\ej\(^g,j^^'6e^'''g^,,, +g,.e>^''e^'''5g,,,> + 7^{x)e^^^5e 

\ I d'-x^\{e'-G^^d-^lG^P5d;^ +e'^G^''5g,.+ ii{x)d^'^5d-l) 



using (53). Using partial integration J d^'^x \/\G\ V = and VG = we obtain 



SS 



-J d'''xy/\G\6A^,{G''^'G''^Vp{e''e;^^) - V,{e-''7]ef"')j + 6<f)' 6^J^, (^y/\G\G^'d^(P 
- j (M,(G'^^G"^^V,(e-^^^,i) - \G\-'^^d,i\G\'^^e-'^r]9P^)') + 5<P' 5,Pg^ 



- / d'''x^y\G\ MJG'"^G"'^Vp(e"^ 



using (58) and (53) in the last steps. This gives precisely the equations of motion (63) and 



(61). Comparing with (56), it may seem puzzling that □gx'^ = must be assumed as well 



in order to get (63). This can be understood as follows (assuming g = G): Since the x"" are 
Cartesian coordinates on M^, Ogx^ = is actually a consequence of Dg^* = 0, expressing 
the fact that Ai C is a minimal surface. Thus from the GR point of view, x'^ is gauge 
fixed. The physical degrees of freedom in are transmuted into 9^'^, which satisfies the 



independent equation of motion (58). 



Finally the propagation of the U{1) degrees of freedom contained in 9'^'^ resp. 6Afj, is 
governed by a Maxwell-like action, which for flat embeddings is precisely the U{1) gauge 
theory discussed above, and for nontrivial backgrounds takes the same formp^ as for the 



nonabelian components discussed in section 3.3 



Remark on extra dimensions. The discussion so far applies to noncommutative spaces 
Ai'^ with arbitrary even dimension d = 2n. Even though the primary interest is on 4- 
dimensional space-time, this is important because any realistic scenario is likely to involve 
extra dimensions, such as Ai'^ = x K"^"^ C ]R^° where K"^"^ is some small compact space. 
We will see that the nonabelian gauge fields obtained on Ai"^ are a priori always su{N)- 
valued, which must be broken spontaneously to obtain some more realistic gauge group. 
A natural way to achieve this is through (fuzzy) extra dimensions, cf . 33 -35 . Such extra 



dimensions may also play an important role in the quantization and will affect the scaling of 
various parameters such as the effective gravity scale. These are topics for further research. 
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this was elaborated by A. Schenkel, unpublished 
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Relation with string theory. The IKKT matrix model was proposed originally as a 
non-perturbative definition of type IIB string theory on M^'^. Our results are indeed very 
reminiscent of D-branes in a -B-field background 36 : Aig could be interpreted as a brane 



in M with open string metric ~ G^j,, while g^i^ corresponds to the closed string metric, 
and is essentially the B field (more precisely B + F, absorbing the U{1) gauge field). 
There are also other solutions of the matrix model consistent with this interpretation, such 
as superpositions of branes. Graviton scattering has been studied from this point of view 
in matrix models, and a relation with supergravity has been conjectured; for an incomplete 
list of references see fTll2ll4l[6l|37H42l and references therein. However, most of this work was 



focused on the interaction of different block-matrix configurations, thereby probing the bulk 
metric and the relation with lOD supergravity. NC branes have been considered only for 
geometries with a high degree of symmetry. The essential new point here is to consider curved 
NC branes in the matrix model and to study the effective gravity on such a 4-dimensional 
brane solution, which does not require any lOD physical bulk. This is indeed possible, and 
it opens the possibility to address 4-dimensional gravity through this simple type of matrix 
model. Thus the strength of string theory (notably the good behavior under quantization) 
seems preserved while the main problems (lack of predictivity) are avoided. In particular, 
the matrix model should be viewed as background-independent, since physical space-time 
emerges dynamically. 



Related aspects and results. The topic of membranes in matrix models has a long 
history, cf. (43-45 . There has been extensive recent work on the noncommutative geometry 



and physics of such "special" NC branes, starting with the fuzzy sphere 32 , which generalizes 
to fuzzy CP"^ 10 and to arbitrary quantized coadjoint orbits, g-deformed versions 46 , etc.. 



Other examples include the fuzzy torus 47 , and various examples with less or no symmetry. 



cf. 48 , 49 and references therein. These and similar spaces may play a physical role as 



compactified extra dimension. Matrix models with additional (quadratic or cubic) terms are 
known to admit such homogeneous spaces as solutions, cf. [7,10 13,39 

" where NC f/(l) 
A relation 



The present framework can be seen as a realization of the ideas in 
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gauge fields were interpreted as gravitons; this will be explained in section 4.2 



between noncommutative gauge theory and geometry resp. emergent gravity has also been 
advocated in |15lfT6ll54 



relating in particular self-dual NC Maxwell theory with self-dual 
gravity. These general ideas are clearly similar in spirit and related to the results presented 
here, although the precise relation is not always clear to the author. Here we restrict ourselves 
to the case of matrix models, which shows the need to consider nontrivially embedded branes 
as well as quantum effects. It is well-known that such matrix models are very rich and admit 
geometrical phenomena (change of dimension, orientifolds, etc.), see e.g. 50,51 . For other 



related work see e.g. |55lp6], and 157] for a different approach to NC surfaces. 
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3.3 SU{n) gauge fields coupled to gravity 

To avoid notational conflicts, we denote the basic matrices with in this section, governed 
by the same matrix model as above 

Sym = - ^Tr[Y'', ,Y''']r]aa'Vw ■ (75) 
However we now consider a matrix background corresponding to n coinciding branes 

= / a = /i = l,2,...,2n, 

\Y' J \(t)'^ln, a = 2n + i, i = l,...,D-2n. ^ ' 

We want to understand general fluctuations around this background. Since the U{1) com- 
ponents describe the geometry, we expect to find SM(n)-valued gauge fields as well as scalar 
fields in the adjoint. It turns out that the following gives a useful parametrization of these 
general fiuctuations: 



wherJ^ 



parametrize the SM(n)-valued gauge fields resp. scalar fields, and A" denotes the generators 



of su{n). This amounts to the leading term in a Seiberg-Witten (SW) map 36 , which relates 
noncommutative and commutative gauge theories with the appropriate gauge transforma- 
tions. This SW parametrization can be characterized by requiring that the noncommutative 
gauge transformation 5„c^'* = ^[1^",A] induces for the su{n) components = A^q,A" and 
the ordinary gauge transformations 

6ciA^ = i[A^,A]su(n) + d^A{x) +0(6) 

5,i¥ = i[^\A]suin) +0(6'), (79) 



cf. 58 . Here the subscript [$*,A]su(„) indicates that only the commutator of the explicit 
su{n) generators is to be taken, but not the 0{6) contributions from the Poisson bracket. 
This strongly suggests that the matrix model action expressed in terms of these should 



reduce to a conventional gauge theory in the semi-classical limit; this was verified in 17 



In our context, the SW map (77) can be understood geometrically in a very simple way. 
Recall that the U (1) sector (i.e. the components proportional to 11„) describes the geometrical 
degrees of freedom: in the geometrical limit, X" = {X'^, 0*) become functions 

= {x^',(f)\x)) (80) 



*We drop sub-leading terms, such as dp^^} for the nonabehan scalars. 
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on A4 which describe the embedding of the 2n- dimensional brane A4 C M^. Then the 
one-form A^dx^ together with the Poisson tensor determines a tangential vector field 

Af.e^' = A^e^'^d^ = A^d^ G TpM ® su{n), (81) 

whose push-forward in the ambient space 

A'd.x" ^ A''{6i:,d,(P') (82) 



coincides with the fluctuations = — X" of the dynamical matrices in (77) (for van- 
ishing $*). This provides the link between gauge fields and "covariant coordinates" X"-. 
The nonabelian scalar fields can be thought of as coordinate functions of Ai embed- 
ded in further extra dimensions corresponding to su{n). They behave as scalar fields, but 
might contribute to the background geometry if they acquire a non-trivial VEV, completely 
analogous to the U{1) components (f)\ 

After this preparation, we can write down the effective action for su{n)-valned gauge 
fields A^ on general A^^" C in the matrix model (75) in the semi-classical limit: 



Sym[A] ~ ^ j d^''x^J]G^\e''G^^''G'"'\i{F^,F,,,.) - Snc 
where = dfj,Ai, — d^A^ + i[A^, A^] is the S'u(n)-valued field strength, and 



(83) 



ri{x)tTFAF. (84) 



n=2 Aq 

~2^ 

The last line holds for 4-dimensional branes, and can be seen using 

^(F A F),,,J^'''9P'' = {F^Jn{Fpjn + 2F^,Fj^''eP'' (85) 



and 6 A 6 = i]e "9 A 9, see 17 . Note that g^y enters the "would-be topological term" S 



NC 



through 7] resp. the antisymmetric matrix 

Qur, ^ G''Pgp^{y)9^"^ = -e" G''P9;^G'''^ = ~9''\ (86) 

This result is non-trivial; for the 4-dimensional matrix model it was first obtained through 
a direct but rather non-transparent computation of the action 17 , which requires the 2nd 
order Seiberg Witten map. For the general case a proof based on the conservation law 
(102) was given [2l!|, which is not only much simpler but also establishes the corresponding 
Yang-Mills equations of motion at the quantum level. 
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It follows that the (bare) 5u{n) Yang-Mills coupling constant is given by 



9ym = ^9'- (87) 

The factor e'^ is reminiscent of a dilaton, however it is not an independent field here but 
determined by the geometry. This implies that e~°" must be nearly constant in any realistic 
solution, which in turn implies that e~'^ ~ A^^^ must be nearly constant. In the case of a 



(nearly-) self-dual symplectic form as discussed in section 4.1, we have rj = e'^ and 

^y^^i^^ = t^J e"tr(ti4x|G^.|i/2G"^M'G'-'F^,F^,,, + 2F A F 

which is essentially a selfdual combination. Recall also that the are su(n) gauge fields, 
hence any realistic gauge theory must descend via some symmetry breaking mechanism such 



as in 34 before phenomenological conclusions can be drawn. 



Nonabelian scalars. It follows easily along the lines of section |3.1| that the geometrical 
action for the su{n) -valued scalars $^ in the matrix model at leading order is 

r , 1 

2(2^0V 

(89) 

where -D^$* = (9^$* + i\A^^ $']. It is worth pointing out that in the case D = 10 and upon 
adding suitable fermions as in the IKKT model [l], this will lead to the analog of = 4 
SYM theory on M^, or on a more general geometry. The SUSY transformations in general 
act non-trivially on the geometry, as indicated below (20) . 

One might question why only the f/(l) components should be considered as geometrical 
degrees of freedom, in contrast to the SU{n) components. The special role of the trace-[/(l) 
components lies in the fact that they couple universally via commutators to all other fields. 
This universal coupling is responsible for gravity, as shown explicitly. In contrast, the SU{n) 
components of nonabelian gauge fields couple primarily through their nonabelian nature. 



3.4 Fermions 



The most obvious (perhaps the only reasonable) action for a spinor which can be written 
down in the matrix model frameworlf^ is 



1 



(27r) 



(90) 



^^In particular, fermions should also be in the adjoint, otherwise they cannot acquire a kinetic term. This 
is not incompatible with particle physics, see e.g. 



34 



35 
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dropping the nonabelian scalars for simplicity (see below). This should be added to the 
bosonic action The Euclidean version involves the obvious replacement ^ — ?• Here 
Ta defines the D-dimensional Clifford algebra. We defined the "local Clifford generators" 



7^(x) = d^x'^Ta, ri^) = e'"/' lu (91) 
which satisfy the Clifford algebra associated with the metrics gfj,u{x) and G'^'^{x) on Ai, 



Thus the (matrix) Dirac operator can be written as 



(92) 



(93) 
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Ip"^ does not quite match with the standard covariant Dirac operator 



where uf'' 



2^ 



V^eJ^ is the usual spin connec- 



As pointed out in 

for spinors Ip^^^J^ ij^e^ [d^ + T.i,^u^^) ^ 
tion, and T^ab = |[7a,7fe]- While the explicit derivative term in (93) is essentially the same, 
the spin connection is trivia^ in the matrix coordinates x^. This is consistent because the 
spinors are those of the ambient space MP . To understand this in more physical terms, re- 
call that the spin connection determines how the spinors are rotated under parallel transport 
along a trajectory. However, can always be eliminated (via a suitable gauge choice) along 



an open trajectory, so that the conventional kinetic term essentially coincides (93) in the 



point-particle limit. Hence the trajectory of a classical fermion with action (90) will follow 
properly the geodesies of the metric G^y, albeit with a non-standard gravitational "spin- 



dragging". This means that holonomies here will be different than in General Relativity, 
and the gravitational spin rotation for a free-falling fermions might provide a characteristic 
signature for the framework of matrix models. For special geometries, ij) coincides with the 
usual Dirac operator on curved manifolds, e.g. in the case of the fuzzy sphere 
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We will see below that upon integrating out the fermions, an induced gravitational action 
is obtained which induces the expected Einstein-Hilbert term J dl^x h? R[G] , albeit with an 
unusual numerical coefficient, as well as additional terms. All this shows that (90) defines a 
reasonable action for fermions in the background defined by G^y. 



Supersymmetry. The presence of fermions in the matrix model is important not only to 
describe physical fermions but also because they allow supersymmetry. The IKKT model [l] 

SjKKT = -^"^^ {^^[X'',X'']{X-' ,X'']r,,a'Vw + ^W^JX^vl/]') (94) 



^^Incidentally, the present Dirac operator on submanifolds is the analog of the operator used in Witten's 
proof of the positive energy theorem j59j . 

^^for massless particles, the geodesies are independent of possible confornial factors. Masses should be 
generated spontaneously, which is not studied here. 
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with D = 10 and \l/ a Majorana-Weyl spinor is singled out by the existence of the following 
maximal supersymmetry 



5(2)^ = ^, 6'^^^X^ = 0, (95) 

where Fab = |[ra,rf,], and e,^ are Grassmann- valued spinors. This model is obtained by 
dimensional reduction of A/" = 1 super- Yang-Mills theory on M^'^ to a point. 

To recover spacetime supersymmetry, we split the matrices again into 4 + 6 dimensions 
X" = Then the 4-dimensional Moyal-Weyl quantum plane M.g is a (BPS) solution 

of the generalized matrix equations of motion, embedded as 

X'' = X'^, /i = 0,...,3, 
(j)' = 0. (96) 

All previous geometrical considerations can be generalized, except that the matrix model 
now contains scalar fields 0*(x). For simplicity we only discuss the case of with flat 
embedding d^cj)' = 0. If we set ^ = ^O'^'T^^e following (l) and recall X^' = - O^'^A^ and 



(91), then the combined transformation S = S^^^ + 5*^^^ takes the form 
6(1)' = ie-f5^i^, 

6 A, = le-^'C^eY^ (97) 

where tf"" = \ [7'', 7''], A^- = ^[Aj, A^], and = 7^, Tg+i = Aj75 where A^ generates the 
5*0(6) Clifford algebra. The constant e°" factors can be absorbed by rescaling the fields. This 
indeed amounts to the A/" = 4 supersymmetry on M^. In the case of general NC backgrounds, 
the SUSY transformation will also act on the metric Gy,i,[x). This can be viewed as a 
supersymmetric form of emergent gravity, which should be worked out elsewhere. 

3.5 Symmetries and conservation laws 

The matrix model is invariant under the gauge symmetry (|4]) as well as a global rotational and 



translational symmetry. The gauge symmetry has been identified geometrically in section 3^ 
as invariance under symplectomorphisms. The SOiD — X^ 1) resp. SO{D) rotation invariance 
is realized as X° — )■ A^ X^ supplemented by suitable action on the spinors, and translations 

cLCtj cLS 

X" ^ X'' + c"!, G M. (98) 

Together these form the inhomogeneous Euclidean group i.e. the Poincare group. As usual, 
these global symmetries should lead to conservation laws. 

Now consider the semi-classical case. Taking advantage of this rotation and translation 
symmetry, one can choose for any given point p G adapted coordinates where the brane 
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is tangential to the plane spanned by the first 2n components, i.e. 9^j0*|p = 0. Then the 
embedding metric satisfies 

= d^5g^,X = 0- (99) 

We denote such coordinates as "Normal Embedding Coordinates" NEC from now on; they are 
preferred matrix coordinates ~ X^. This simplifies the analysis of general branes, allowing 
to reduce many considerations to the case of trivially embedded branes. For example, it is 



now easy to see that the matrix model action for the U{1) fields can be written as in (66) in 
terms of 

47]{x) = {x", x''}{x'\x''}r]aa'Vbb' = G"^" (x) (7^, (x) . (100) 

This is so because ri{x) can be viewed as a SO{D) scalar, which in normal embedding 
coordinates reduces to e" G^^{x)iri^^. The covariant version is as above. 

If the effective metric coincides with the embedding metric G^i, = g^^ as discussed below, 
then these NEC are "free-falling" Riemannian normal coordinates in GR, i.e. G^y\p = r]^,^ 
and (r|^^)|p = 0. Local Lorentz invariance is then manifest in these coordinates. However, 
an important distinction with GR becomes obvious: the extrinsic curvature of the brane 
Ai C does play a role in the matrix model (at least in the "harmonic branch" discussed 
below), due to the brane tension. For example, even though a cylinder embedded in is 
intrinsically fiat, it is not a minimal surface. 



3.5.1 Tangential conservation law 



Motivated by the translational symmetry, one can consider the following "x-dependent tan- 
gential" variation 

X" ^X'' + 5X\ (5X" = {X\ [X^ eb(X)] } (101) 

which defines a measure-preserving transformation on the space of matrices. The corre- 
sponding Dyson- Schwinger equations have the form of a conservation law 



[X\T'^'']7^aa' = Q. 



(102) 



This provides useful non-perturbative information which is not restricted e.g. to BPS sectors. 
The bosonic action ([T]) leads to the matrix - "energy-momentum tensor" 18,61-63 



It is easy to verify this using the matrix e.o.m., but (102) is more useful because it involves 
only tangential tensors on M.. Indeed the semi-classical geometrical limit of T"** is given by 



T 



ab 



-2G'^'^(x)5^x"9^x' + 2'q''^7]{x) 



(104) 



which for g^y = G^i, is essentially the projector on the normal bundle of C M.^. As 
elaborated in 21 , the conservation law (102) then leads to the following equivalent equations 



(105) 
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(106) 
(107) 



where 6'^^ is anti-symmetric and defined in (86). The last form coincides precisely with 



the e.o.m (63), which therefore can be trusted at the quantum level, and it holds more 
generally e.g. in the presence of a mass term or a more general potential. Note that the 



D relations (102) reduce to 2n tangential ones, while the "transversal" components give no 



new information. 

This matrix conservation law also applies to the nonabelian SU{n) components. Then 
(102) gives precisely the Yang-Mills equations of motion for SU{n) gauge fields coupled 
to gravity obtained from (83). This seems to be the best way to establish the action for 
nonabelian gauge fields I21j. 



4 Special geometries and perturbations 

The most obvious solution of the above covariant equations of motion is of course M^", 
corresponding to the semi-classical limit x°' = {x'^, 0) of the Moyal-Weyl quantum plane M^". 
Its effective metric G^'^ is flat. In the supersymmetric case, these are BPS solutions |l] and 
therefore protected from quantum corrections. To understand other, less trivial solutions, it 
is illuminating to introduce the following tensor 

Then the effective metric can be written as 

which is the reason for choosing the above normalization. J'^ satisfies the following properties 

trj^ = -ie-^r] (110) 



due to the anti-symmetry of 6^^. In 2n = 4 dimensions, we have the additional relations 
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detj = 1 

J'' + J~^ = ItijH (111) 



using (52). We can now extend the discussion in section 2.2.1 on (anti) self-dual symplectic 
structures and Euclidean versus Minkowski signatures to the general case. 
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4.1 Self-duality and G^i, = g^^. 

We focus on the 2n = 4 - dimensional case. Since everything is now formulated in tensorial 
language, the same arguments of section 2.2.1 apply at any given point p E Ai. Assume 
first that Ai is Euclidean. Then we can diagonalize the embedding metric at that point 
g/iulp = Sfj,iy, and bring the Poisson tensor into canonical form 



9 (a dx dx ± a dx dx 



(37) at p E Ai using a suitable 50(4) rotation. This leads tc^* 

G^"" = diag(a^ 



— 2 2\ 



at p E M. 



(112) 



(113) 



Noting that •k^dx'^dx^) = dx^dx^ where -k denotes the Hodge star defined by et^'^p'^ and g^^i^ 
on Ai'^, it follows again that the corresponding symplectic form is (anti-) self-dual ((A)SD) 
if and only if 



•k CO 



±L0 ^ 



a 



1 resp. e "7] = ! ^ G^^ = g^iu ^ 



-1. 



These are tensorial statements. The second statement follows from 

1 



V 



(a^ + a^^) > 1- 



:il4) 



(115) 
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Note that JT"^ = — 1 means that JT" defines an almost-complex structure, and Ai becomes 
an almost-Kahler manifold. It is easy to check that then (63) reduces to V^^^J = 0, which 
follows from -ku = ±u resp. d-kuj = ±du = 0. 

In the Minkowski case, completely analogous statements apply once we adopt complexi- 
fied Poisson structures with imaginary time components, having in mind a Wick rotation 

— )■ it as discussed in section 2.2.1 If we accordingly define e^^^^ = then 

-k (jj = <^=^ a = 1 resp. e~'^?7 = 1 <^=^ G 

where 



= -1 



(116) 



9 ^{iadx^dx^±a ^dx^dx'^). 



(117) 



Notice that = —1 as in the Euclidean case determines the reality structure of 9^'^ , and is 
replaced by the quartic equation (111) in the general case g^^, ^ G^y. This ensures that the 
physical metric 

G'^'^ = diag(-a^ a-^ a-^ a^) at pEM (118) 

is always real. 

It is important to note that self-dual Poisson structures always solve the equation of 
motion (63), which reduce to V^9~^ = for e°" = rj. Such (A)SD closed 2-forms 9~^ with 



^^Of course this holds only at p and does not extend into a local neighborhood of p e A^. 
^^A reader not comfortable with such a complexified 9'^'^ may prefer to stick with the Euclidean case and 
postpone the issue of Wick rotation. 
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constant asymptotics 6~^{x) — )■ 6^^ as a; — > oo always exist on (suitable) asymptotically 
flat spaces. This can be seen by interpreting as sourceless electromagnetic field with 
constant field strength at infinity: we only have to solve rf^F = 0, F = dA with constant 
asymptotics F — )■ F as r — )■ oo, and define to be the (A)SD component of F . Therefore 



G 



fii/ 



0, 



□c 







(119) 



provides a particularly interesting and transparent class of solutions without matter. The 



bare action (66) then has the form of a brane tension using (114) 



S 



A'' 

%2 



9^ 



(120) 



This has the same form as the one-loop vacuum energy, and the two may approximately 
cancel provided Aq coincides with some effective cutoff. In that case, the total effective 
brane tension may indeed be small. However, even though minimal surfaces are in general 
not BPS, they are protected and in fact stabilized by quantum corrections. 

In the remainder of this review, we will mostly assume that uj is (anti-) selfdual to a 
very good approximation, i.e. G^j^u = g^^u- Most of the known solutions of matrix models 
(possibly with additional quadratic and cubic terms) such as fuzzy CP^ etc. indeed satisfy 
this condition. A more in-depth study of this class of geometries is given in 79 . However, 
we now study small perturbations of this (A)SD case. 



4.2 Metric perturbations from U{1) gauge fields 

Using the above results, the most general geometry of the matrix model can be viewed as 
perturbation of an (A)SD background G^j, = g^^ through "would-be" U{1) gauge fields 
with field strength F^^, leading to a metric perturbation 

G^u = G^y + h^y (121) 

elaborated below. These U{1) gauge fields result from tangential variations of the matrices, 
while transversal fiuctuations of the matrices lead to variations of the embedding metric g^j,^. 
This appears to be the most useful parametrization. 

Hence consider some M. C M"^ with embedding normal coordinates at p G such that 
Qnu = diag(±l, 1, 1, 1) at p e A^. We can assume that the Poisson tensor 9^^^ has the form 



(117) resp. (112), leading to an effective metric 



G^y = diag(±a^ a'^, a'^, a^) (122) 



at p G Ai as in (44), which in the case of (A)SD Poisson structure satisfies G^^ = g^^. Now 
consider small fiuctuations of the tangential matrices around such a configuration. They 
can be parametrized as 

X^" = X''- 9'"'Ay{x) , (123) 
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and the A^, were interpreted in section 2.2 as U{1) gauge fields. On the other hand, we saw 
in section |3.2| that they can be interpreted as perturbation of the effective Poisson tensor 

(124) 



(pl on M, 



0-J + 



flU 



where = d^A„ — d^A„. This in turn leads to a perturbation of the metric (121) on M., 
with 

V = ^^((:f-'):;>M'/. + (^"')^^^v-^((^-')^Fp.G'^'^)G^.) (125) 
where p = lOZul^^"^ and 



1 











—a 












—a 












a 























I 



(126) 



To be specific we assume the Euclidean case here, where 9^^ has the canonical form (44). 
As a check, observe that h = G^'^hni, = consistent with (54). This interpretation of 
the noncommutative U{1) gauge fields in terms of a metric perturbation was first given by 
Rivelles 
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He observed that if the gauge fields solve the Maxwell equations, then the h 
are actually gravitational waves resp. gravitons. To see this, consider the case of Ai 
so that the linearized Ricci tensor is given by 



4 



2 ^ ^ 2 



5 '''/i-y ■ 



(127) 



Now recall the equations of motion (65) (equivalent to the Maxwell equations for F^j,) in the 
preferred (matrix) coordinates, which together with (59) gives 







For the metric fluctuations this gives 



Gna \ G 



1 



- = 0, 

and the expression for the linearized Ricci tensor simplifies as 



S'l'ii'y 



Now the Maxwell equations for the U{1) fluctuations on 



imply as usual 



fii/ 







0. 



(128) 
(129) 

(130) 

(131) 
(132) 
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Together with (130) and (125), it follows that these U{1) metric fluctuations are Ricci- 
flat [14j, 

R^.u[%. + h^u] = 0. (133) 



This can also be seen as a consequence of Oqx^ = 0. Once we establish that the corre- 
sponding Riemann tensor is indeed non-trivial, these are genuine gravitational wave; 
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m 



harmonic gauge (129)), parametrized through on-shell degrees of freedom of U{1) gauge 
fields. This provides an interesting relation with GR. 

To make this more explicit, consider the case G = g i.e. a = 1. Using the usual 
parametrization of F^,^ in terms of Ei and B^, the corresponding metric perturbation is 



h 



/ii/ 



[ (Bs-Es) 

{B2 - E2) 

-{Bi-E^) 
\ 



{B2 - E2) 
-{B, - E,) 


(Bi - EA 



-{B,-E^) \ 

(5i - E^) 

-{B3 — E3) {B2 — E2) 

[B2 — E2) (-B3 — E^) J 



(134) 



This vanishes for self-dual fields Bi = Ei as expected due to (114) (only for a = 1!), and 
is non-trivial for the anti-self dual fields as we verify below. In the Minkowski case, the 
corresponding result for 



e^'^ = e 



/ 



1 

\ i 





-1 








/ 



and F, 



flU 



\ 



IS 



h 



\ 



-{E3 + B;) 

-l{B2 + E2) 

i{B^ + E^) 




-l{B2 + E2) 

-{B-s + E,) 


Bi + E, 



iEi 


iE2 iE-i 


\ 


-iEi 


B3 —B2 




-iE2 -B3 


Bi 




-iEs B2 


-Bi 


) 


t{Bi + El) 


\ 







Bi + Ei 




-(^3 + ^3) 


B2 + E2 




B2 + E2 


Es + B; J 





(135) 



(136) 



which vanishes foi^ -kF = F. The imaginary hoi are misleading, since the physical metric 
in real coordinates is manifestly real (118). The above form is written in "partially Wick 
rotated" coordinates, and the real physical form is seen in coordinates where hoi = 0. This 



is obtained after a transformation — )■ + 



x] 



X 



+ ite 



x] 



which amounts to the 



required local Lorentz boost to achieve the normal form (117). In case of doubt, one can 
first work in the Euclidean framework and do a formal Wick rotation in the end. 

Hence on a self-dual background ^f^y = G^^, only the anti-self dual perturbations F^j, de- 
scribes non-trivial gravitational waves, while the self-dual perturbations do not lead to any 
metric deformations at least at the linearized level. However they are non-trivial perturba- 
tions of 6^'^. Hence they seem to be completely sterile fields which decouple from everything 
else. This suggests that they should be integrated over. That would in particular suppress 



Lorentz- violating effects due to 9'^", such as in the term J d'^x R^upa 



e^^QP'' found in 



On 



It was argued in 14 that these metric perturbations generally correspond to pp waves. 



^^recall that we use e^^"'^^^ = i in the Minkowski case. 



30 



the other hand if 6'^'^ is taken to be real in the Minkowski case, then one finds 2 independent 
physical gravitons as in GR. 

In summary, configurations with g^^, = G^y can provide the required geometrical degrees 
of freedom for full-fledged gravity assuming D = 10 ^O], while perturbations due to the 
"would-be" U{1) gauge fields are ac?(iztzona/ graviton-like degrees of freedom on such a back- 
ground. On the other hand they can also be interpreted in terms of the scalar fields (j{x) 
and ri{x). Clearly more work is required in order to disentangle and fully understand the 
significance of these degrees of freedom. 



Plane waves and linearized Riemann tensor. We verify here that the Riemann tensor 
of the above metric perturbations is non-trivial. Consider in particular plane waves 

A^{x) = A^e'^- (137) 

such that F^u = i{kfj_Ai, — k^A^) ^ 0. Then 

V = ie{k^A, + KA^)-id{k^A, + kyA^)+iG^,{eP''kpA^) (138) 

where 

h = {J-XK- (139) 
Observe that the term i{k^Ay + kjjA^) has the form of a diffeomorphism generated by Ay, and 
therefore does not contribute to the Riemann tensor. Using (138) for the metric fiuctuations 
gives 

-^o-/]iv7 ~ ~ 2^ (^i^fi^cT ~ kf^kcr){kyA^ — Ayk^) -\- (^k^A^^ — k„An)(^kykj — k^k^)^ 

2 i^i^^i^^^cri kykfjG k^k^G ~\~ k^kfjG ^y)(d^ kpA^-) ^ 0. (140) 

Thus the U{1) metric fiuctuations correspond to propagating gravitational waves on with 
2 independent degrees of freedom, i.e. physical gravitons in harmonic gauge 
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Linearized coupling to matter. Here we briefiy discuss the coupling of these "would- 
be" f/(l) de grees of freedom to matter, corresponding to perturbations of Gf^y — 9fj,u- 

The 

coupling of hfj,y = G^y — g^y = SG^y is given as usual by 

55^attcr ^SnJ d^xVGSG^yT^'' (141) 

where T^y is the energy- momentum tensor of matter. We will focus on the case of static 
matter distributions with Tqo = p(a;), Tjj = Tqi = 0. Then the linearized metric perturbation 
(121) couples as 

^•S'matter ^ Snj d'^x h^yT^'" = 8n j d'^xhoop{x) 

d^x {E-mE + B- mB) (142) 



^the restriction to (A)SD field strength is not seen here because there are no (A)SD plane waves. 
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in the Minkowski case, replacing {—E^ — B^) — )■ [B^ — E-^) in the Euchdean case. Hence from 
the point of view of electrodynamics, matter behaves like a dipol^^ with dipole strength 

rfiE = -8710(0, 0, p(x)) = ms, (143) 

pointing in the direction determined by The total dipole moment is given by the mass 
M, 

Mj^ = -Sit M 063 = Mb, (144) 

Note again that only the anti-selfdual components of the would-be U{1) gauge field couple 
to matter. In particular, observe that the total effective dipole moment of e.g. a system 
of stars rotating each other is time-independent and pointing in a constant direction, hence 
they can radiate only quadrupole and higher multipole modes in /i^^. 

Hence matter leads to a small localized dipole-like deviation from 6'^'^{x) being self-dual, 
leading to g^j,^ ^ G^y near matter. However, this will be suppressed by powers of ^ resp. 
-f^i and is therefore presumably negligible for most purposes. We will therefore ignore the 
would-be U{1) gauge fields, and focus on geometries with G^^ = g^y. 



5 Quantization 

Perhaps the main motivation for the matrix model approach is that they provide a natural 
non-perturbative concept of quantization, defined in terms of an integral over all matrices 

ciX"(i^e"^I^'*l (145) 

and similarly for correlation functions. In particular, the IKKT model thereby promises to 
provide an accessible quantum theory of gravity (as well as the other fundamental interac- 
tions), because it can be viewed equivalently as A/" = 4 NC SYM on Mg. This allows to 
circumvent many of the technical problems in quantizing GR. In the Euclidean case, it was 
shown that this type of integrals exists for finite (in D > 3 dimensions) (67,68 , even 



for purely bosonic models; see also 69 for some non-perturbative results. Nevertheless, in 
order to obtain a well-behaved quantum theory, it is probably necessary to restrict to the 
IKKT model with A/" = 4 supersymmetry resp. D = 10, or closely related models; this 
will be explained below. Remarkably, D = 10 happens to provide precisely the degrees of 
freedom needed to realize the most general 4-dimensional geometries through embeddings 
A^^ C ]R^° |70|. Even though the mechanism is quite different, it is worth recalling that 



single-matrix models have been used as a non-perturbative definition of 2D gravity 71 . 

In the following we study the quantization around some NC background as discussed 
above, which amounts to a "condensation" of the matrices. This is a strong assumption, 
which would imply that the f/(oo) gauge symmetry is spontaneously broker However, 



^''This is strongly reminiscent of the picture of [65[|66| where open strings act as dipoles in the presence of 
a -B-field background. 

^^gauge fields can thus be viewed as analogs of Nambu-Goldstone bosons, where the gauge symmetry is 
realized non-linearly. This is a beautiful insight provided by NC gauge theory. 



32 



many of our results should hold also under somewhat weaker assumptions, notably {X^) = 
{[X^^X"]) = but {[X^',X•'][X'',XP]) ^ 0; more specifically, O^"" might not have long-range 
order but short-range fluctuations and should be averaged. This may also be important for 
suppressing Lorentz- violating effects. 



5.1 UV/IR mixing in noncommutative gauge theory 



section 



2.2 



Consider first the matrix model from the point of view of NC gauge theory on as in 
A general matrix resp. field on can be expanded in a basis of plane waves. 



<P{X) = J ^ 0fc e^'^"^' e Mat(oo, C) - (146) 

where 0^ G C is now an ordinary function of A; G M^. The free (quadratic) part of the 
action is then independent of 6^^", but the interaction vertices acquire a nontrivial phase 



factor e2 ^'<j where fc^ denotes the incoming momenta. The matrix integral (145) 
then becomes an ordinary integral J dX"- = jHrf^fc, which can be evaluated perturbatively 
in terms of Gaussian integrals, similar to ordinary QFT. This leads to the well-known Filk- 
type Feynman rules 72 for field theory on Mg, where planar diagrams coincide with their 



undeformed counterparts, while the non-planar diagrams involve oscillatory factors. All 
this can be made rigorous on fuzzy spaces, where the matrices are finite-dimensional. The 
reason for recalling these steps is to emphasize that these rules are consequences of the basic 



definition (145), there really is no choice in the framework of matrix-models. 

In particular, the loop integrals are generally divergent in spite of the existence of a 
fundamental length scale A^c- The reason is that area rather than length is quantized on 
M.g. This implies that a UV-divergence in some direction fc^ ^ A^vc is necessarily associated 
with an IR effect in another direction; this is the essence of UV/IR mixing [73|. Assuming 
some cutoff A in the loops, one finds new terms in the one-loop effective action for the trace- 
U{1) components, which are divergent for external momentum p — )■ 0. In the example of 
scalar fields coupled to external U{1), they have the form 

+ i(9F(p))(eF(-p))(A:, - ip.pAjj + (47 - 301og(^)))) (147) 



where 



1/A 



^^The generalization to the Minkowski case is mathematically more subtle, but formally the same derivation 
applies and preserves the basic oscillating character of the vertices. This should be well-defined at least in 
the case of (softly broken) J\f = A SYM, consistent with the concept of Wick rotation discussed previously. 
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The first term amounts to a renormahzation of the couphng, but the second term is unusual 
and singular as p — )■ 0. These IR divergences become worse in higher loops, and the models 
are probably pathological as they stand. There is a notable exception given by the A/" = 4 
SYM theory on M^, which is nothing but the IKKT matrix model for D = 10. This is 



(almost ) the unique model which has no UV/IR mixing at one loop, and is arguably well- 



defined and finite to any order in perturbation theory just like its commutative cousin |24). 

Now focus on the A/" = 4 model with soft SUSY breaking terms, leading to a physical 
cutoff A = A4 given by the scale of = 4 SUSY breaking. Then the above "unusual" 
term in r$ still leads to physically unacceptable effects from the point of view of U{1) gauge 
theory (e.g. polarization-dependent dispersion relations for the photon [74]), however they 
are perfectly sensible from the point of view of emergent gravity: they amount to an induced 
gravitational action, as explained below. 



5.2 Induced gravity 

We now explain how this phenomenon of UV/IR mixing in NC gauge theory can be under- 
stood in terms of induced gravity. The "strange" new term in the effective action are simply 
induced gravitational termt^^ such as an Einstein-Hilbert term 



Induced gravitational action due to scalars Consider the quantization of a scalar field 



if coupled to the matrix model as in (48). Upon integration, this leads to an effective action 

1. 



e 



dip e 



-5M 



where F,, 



-TrlogAc. 



A standard argument using the heat kernel expansion of Aq gives 77 

1 1 



167r2 



rf^x ( -2A^ - -R[G] A' + 0(log A) 



(149) 



(150) 



which is also the one-loop induced gravitational action due to a scalar field in the matrix 



model on a generally embedded brane /A C 



(147) using the identification (125), in the semi-classical limit 



19 . This is essentially the mechanism 



of induced gravity 78 . It was shown in 19 that this reproduces the UV/IR mixing terms 



pA < Ale 



(151) 



where p is the external momentum. This means that the phase factors e^'^'^^"^^" in the non- 
planar loop integral due to the non-commutativity are small and thus well reproduced by 
the semi-classical Poisson structure incorporated in Gfj_y. 
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which however has the same field content 



^^there is also a particular M — 1 model with that property 
and is closely related to the IKKT model. More general soft SUSY breaking terms (i.e. quadratic and cubic 
terms in the matrix model) should also be admissible. 

it seems likely that the interpretation of UV/IR mixing in terms of closed string modes put forward 
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75 76 is another way of looking at the same thing, since these closed-string modes are related to the 



embedding metric of the brane. 
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Induced gravitational action due to fermions The one-loop effective action due to 
integrating out a fermion in the matrix model is more complicated. Due to the non-standard 
spin connection, one cannot simply use the standard results for the usual Dirac operator. 



The induced gravitational action was determined in 64 using the general formulae for the 
Seeley-de Witt coefficients for fermions. This is still too complicated in the most general 



case; for a Dirac fermion on generally embedded branes Ai^ C with G 
induced gravitational action can be written as 



the 



^ I d'x^\(2A' + ( - ^-R[g] + \d,ad^a 
+le--R[gUpJ'"'9''' + ^(□^x'^) (□3X^)7/,,) + O(logA)). 



(152) 



Here k is the number of components of D-dimensional Dirac spinors. This formula generalizes 
to the case of 2n-dimensional branes. The term {ngX°'){\I\gX^)riab vanishes for on-shell 
(tree-level) vacuum geometries, but may be non-trivial in general, reflecting the embedding 
of C 



3D 



with G^,^g^, 20,64 



A similar result was also obtained for such on-shell (tree level) geometries 
Again, it was verified that this reproduces the UV/IR mixing terms 



the fermionic loops in the semi-classical limit (|151|). There is indeed a corresponding matrix 
model action 



79 , as given below. 



Induced gravitational action due to gauge fields The induced gravitational action 
due to integrating out gauge fields has not yet been computed directly for a general back- 
ground. However, it can be obtained indirectly based on supersymmetry: note that the 
fermionic contribution to the one-loop effective action does not quite cancel the scalar con- 
tribution. Some UV/IR mixing may remain even in supersymmetric gauge theory on M^, 
and the UV/IR mixing is absent only in the case of A/" = 4 supersymmetry 24,25 . This 



strongly suggests that the model should be finite above some scale A to obtain a well-defined 
quantum theory. This is realized by the IKKT model on a NC background, where A is the 
scale of A/" = 4 SUSY breaking above which no divergent terms are induced. 

Hence consider a matrix model with ns scalar fields and nij, Dirac fermions. The Af = 4 
model has ns = 6 and n^, = 2, in addition to the U{1) gauge field A^. Because this model 
is free of UV/IR mixing , it follows that 



-2r, 



6r. 



(153) 



This matches the vacuum energy contribution of 2 scalars (150), but has a distinct Einstein- 



Hilbert term^^ as well as a R66 due to the would-be topological term (85). Note that (153) 



applies to both the SU{n) as well as the f/(l) contributions to induced gravity, because all 
fields in the loop (abelian or nonabelian) couple in the same way to external U{1) legs. 

Thus there will be no induced gravitational action above the scale A of A/" = 4 SUSY 
breaking. What happens below that scale depends on the specific breaking mechanism, 

^^This term can be expressed in terms of intrinsic quantities including 0^''. 

This is the reason why even M = 2 NC SUSY gauge theories do have UV/IR mixing. 
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which could be due to some Higgs effect, fermion condensation, or some other unknown 
mechanism. The R99 may be problematic due to the explicit presence of 9^'^ which breaks 
Lorentz invariance. This term might disappear due to fluctuations of 9^^ , or due to suitable 
SSB such that the contributions from gauge fields and fermions cancel, or by counterterms 



in the matrix model (155). Then the effective gravitational action can take the desired 
Einstein-Hilbert form 

Tgrav = j d'x {K\ + R[G] Kl + 0(log A)) (154) 

plus possibly some additional terms as above. Higher-order contributions in the curvature 
will be suppressed by powers of ^; note that in induced gravity there is no coupling constant 
of negative mass dimension unlike in GR. Moreover, there are no unstable modes in the 
matrix model unlike in GR. Furthermore the strong RG running of the U{1) sector is cut off 
at the NC scale. All this suggests that the quantization of gravity in this model should be 
well-defined and under control. This would be jeopardized in the presence of higher-order 



terms such as (155), requiring and a more sophisticated RG-type analysis. In any case, the 



model certainly provides an exciting new approach to quantum gravity. 
5.2.1 Higher-order terms in the matrix model 



As shown recently 79 , the above induced gravitational actions can be identified with higher- 



order terms in the matrix model. In particular, 

Tr[[X^X^], [X,,X%Xa,X,] ~ j d'x^ -2R + e'^d^ad.a^ 

Tr([X",T^^][X„,T,J+T'^''nX,nX5) ~ j d'x^ {{D - 4)e"nge'^ - 2e^''R) (155) 

assuming (yf^j, = G^^, where OX"' = [X'^,[Xc, X"-]], and matrix indices are lowered with 
r]ab- A matrix expression for e^DgC'' term can also be found. This could be used to cancel 
unwanted contributions in the induced action, and to study systematically the quantization 
of gravity through purely bosonic matrix models. 



6 Macroscopic description: Harmonic branch and Ein- 
stein branch 

In this final section, we want to discuss the low-energy effective gravity theory arising from 
the D = 10 matrix model, trying to minimize any explicit reference to the detailed noncom- 
mutative origin. We will derive equations of motion for gravity coupled to matter in the 
semi-classical limit, assuming that g^^y = G^iy with appropriate (anti-) self-dual 9^^". This 
should be a very good approximation even in the presence of matter. Our starting point 
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is the semi-classical effective action of the matrix model (120) together with the action for 



matter, as well as the Einstein-Hilbert action R[g] induced at one-loop 

S = j d'x^l (AlR - 2A^) + S^,,ter; (156) 
possible additional terms are ignored for simplicity. Now recall 



6 j ^\R = -J ^MQ^'^g,^ 

^^matter = StT / v^T^^^^, (157) 



where Q^" = R^^ — \g^^R is the Einstein tensor, and T^^^ is the energy-momentum tensor for 
matter (recall that matter and fields couple to the effective metric essentially in the standard 
way). The crucial point is now that the fundamental geometrical degrees of freedom are not 
the (yf^i/, but the embedding fields 0* as well as resp. 0^'^ . The most general variation can 
thus be decomposed into variations of g^^ = r^^,^ + 5^0*9,^0' and variations of 9'^'^. The e.o.m. 
for 6'^'^ are satisfied for gfj,^, = G^y to an excellent approximation as explained above. The 
variation of S with respect to the fundamental fields (jf can be written in matrix coordinates 
= X^^ using 5g^y = d^(j)^dy5(j)^ + d^dcjfdijcj)^ as 

^S = j d^x^/\i\5g^yW = -2 J (fx5(j)'d^{^/W\W)dy(t)' (158) 

up to boundary terms, where 

W = SttT^^ - Alg^"" - A^g^"". (159) 
This leads to the equations of motion for 0* 

d^{^\Wd,(l)') = 0, (160) 
which using the identity V^V^ = -^df^{y/gV^) can be written as 

Atng(j)' = {SttT'^'' -Alg^')V^d,(j)' + 87r{V^T^''')d,(f)' (161) 
recalling V^Q'^'^ = 0. This equation has 2 types ("branches") of solutions: 

1. "Harmonic branch" : 

The first class of solutions 

d,{^M'H'"'du<P') = 0, n'''' ^ (162) 
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can be interpreted as solutions of a deformed Laplacian. The prototype of such a 
solution is flat Minkowski space 0* = 0. Remarkably, this is a solution for arbitrarily 
large vacuum energy resp. Ai > 0: the term J d'^Xy/g is interpreted as brane tension 
rather than as cosmological constant. More generally if the vacuum energy dominates 



the matter density, (162) reduces to Og(f)^ fa 0, leading essentially to minimal surfaces. 



As explained in more detail below, this leads to the near-realistic cosmological solutions 



of FRW type 80 , which are stable and largely insensitive to the detailed matter 
content. Remarkably, Newtonian gravity is also obtained due to the brane tension, 
without even using the Einstein-Hilbert action. In that scenario it is natural to assume 
that 

A4 < Apianck (163) 



so that the induced gravitational action and in particular the R99 term (152) are only 
small corrections and pose no strong constraints. Hence the harmonic branch may 
provide a solution of the notorious cosmological constant problem, which is a strong 
motivation for more detailed studies. 

Another attractive feature of this harmonic branch is that its quantization should be 
comparably simple, as the embedding fields 0* are governed by a simple action with 
positive excitation spectrum. On the other hand, it is questionable if the solar system 
precision tests can be met, and more thorough studies are required. 

2. "Einstein branch" : 



The above deformed harmonic equation (162) becomes void if Uni, = 0; then Ai is in 



some sense a tensionless brane. This amounts precisely to the Einstein equations 

Ag^"" + g^"' = 87iGT'"', (164) 

upon identifying 

^=Al § = At (165) 
In particular, the Planck scale should be identified as 

A4 ~ Apianck, (166) 

and the term J d'^x^/g Af appears to be a cosmological constant. Indeed using general 
embedding theorems |70], any solution of the E-H equations can be realized (locally) 
using embeddings in D > 10. This provides a realization of (induced) Einstein gravity 
through matrix models. Since the quantization is defined through the matrix model in 
a background-independent way, certain problems of quantizing gravity such as unsta- 
ble modes are not expected to occur. This should provide a quantum theory of gravity 
very similar to general relativity, along with all the other fields in the IKKT model. 
Singularities should be resolved due to the quantum structure of space-time. On the 
other hand, the usual cosmological constant and its fine-tuning problems appears to 
arise again. However, the different parametrization of the geometry in terms of em- 
bedding rather than a fundamental metric might come to rescue: In the fiat case, the 
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first variation due to transversal 50* drops out (!), and the tangential variations 
are transmuted into variations of tlie Poisson tensor wliicli also drop out as explained 
before. The same argument "almost" applies to near-fiat geometries, which should 
thus be at least meta-stable even in the presence of large vacuum energy. It remains 
to be seen whether the c.c. problem can be solved in this way. 

Hybrid solutions also conceivable where the rank of H-p^u is not maximal. This is expected 
in the presence of extra dimensions M. = x K, where could be a solution of the 
Einstein equations while K could be e.g. a sphere (stabilized by a flux). Indeed the IKKT 
model provides just enough degrees of freedom for Einstein gravity on a 4-dimensional brane, 
but not for higher-dimensional ones. In the presence of additional terms in the matrix model 
such as Tr{X"'X^X'^eabc) or Tr{X°'X^riab), new types of solutions arise. Cubic terms can be 
interpreted in terms of fluxes, and preserve the translational symmetry of the matrix model. 
The prototypical example of a solution in presence of a cubic term is the fuzzy sphere Sj^, 
which is also a solution in the presence of a quadratic "mass" term. 

The solar system precision tests clearly support the Einstein branch, and its realization 
in the present framework should be studied in detail; e.g. the Schwarzschild geometry will 
be discussed in jSlj. On the other hand, the harmonic branch is very interesting because 
it leads to a stable cosmology which is remarkably close to what we see, and to interesting 
long-distance modifications of gravity. We will therefore discuss the harmonic branch in the 
remainder of this paper. Although at this early stage of investigation it does not appear to 
be fully realistic, some of the result might carry over into a modified Einstein branch. 

7 Harmonic branch: gravity from branes with tension 

Based on the general results above, we can start to explore the physical properties of gravity 
in this matrix model. In this section some preliminary and particular results for the "har- 
monic branch" of the model will be discussed, which may or may not be physically relevant. 
These solutions can be understood largely in terms of embedded branes Ai C M^^ with 
tension coupled to matter. We will identify an intuitive and physically appealing mechanism 
for (Newtonian, at least) gravity on the brane, which applies in a more general context. 

7.1 Cosmological solutions 

Consider first cosmological solutions of Friedmann-Robertson- Walker (FRW) type. These 
are characterized by spatially homogeneous and isotropic metrics 

ds^ = -dt^ + aUfdJ:^ dT? = W^rfr^ + r^rffi^ (167) 

1 — K^r^ 

where (iE^ is the metric on a 3-dimensional space with constant curvature, and a{t) is the 
cosmological "scale" parameter. There are 3 possibilities: the space-like 3-manifold is fiat 
(denoted by A; = 0), a 3 sphere 5*^ with k = +1, or hyperbolic space where k = —1. In 



39 



GR, such FRW solutions are obtained in the presence of a homogeneous matter distribution, 
which leads to the Friedmann equations for a{t). 



In the present framework, we will look for solutions describing an embedding Ai C 



)10 



with metric G^j^ = of FRW type. We know from section 4.1 that this is a vacuum solution 



of the matrix model provided /A C M is a minimal surface i.e. Dx" = 0, equipped with 
a (anti-) self-dual nondegenerate 2-form 6~^. The presence of matter will in general modify 
these solutions, however they are expected to be a good approximation provided the vacuum 
energy dominates the energy density p due to matter. Thus we will assume that Ai 
is at least of order eV, say, and increasing Ai will only stabilize the solutions. This is in 
striking contrast with GR, where the vacuum energy resp. cosmological constant must be 
very precisely tuned to Ai ^ 2meV according to the ACDM model. This leads naturally to 
very reasonable solutions which are in remarkably good agreement with observations (leaving 
aside aspects related to the early universe such as the CMB, which need to be re-analyzed 
in the present framework). 

Generalizing the one-parameter solutions in 82 , the following embeddings of FRW ge- 
ometry through minimal surfaces A^^ 



C 



( sinh(x) sin6' COSV9 \ \ 
sinh(x) sin ^ sin 
sinh(x) cos B 

\ cosh(x) 




3 10 



(168) 



where 

and riah = diag(+, ...+, 



TZ{t) = a{t) 



cosipit) 



sin ^{j{t) 

-,+,+). This leads to a FRW geometry with k 
-dt^ + a{tfdY?, dY? = dx^ + smh\x)d^'^ 



(169) 



(170) 



where r = sinh(x). This is harmonic Dx" = provided 

= □(7^(t)sinh(x)cose) 
= Dxc. 

This leads to 



(171) 



3- (a — 1) + d — ip a 
a 



bipd + ipa 

1 . . 

3 QjOc Q ~\~ Xq 





0. 



(172) 
(173) 
(174) 
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Similar solutions for k = +1 were also found which are however unrealistic, while for k = 



only the one-parameter solution in 82 is known. The equations (174) can be integrated as 
follows: 



Da' 



1,2 -2 

a 

d^Xr 



b 
d 



= const 
const > 
■ const. 



(175) 



leading to a 2-parameter family of solutions; the case 6 = was obtained before in [82]. For 
the Hubble parameter this implies 



d^ 



a 
a 



-3d^a-^ + Ab'^a 



-10 



(176) 
(177) 



For the early universe i.e. small a and 6 7^ 0, this leads to a big bounce with minimal size 
ao ~ d^^. There is a transition determined by ci to a coasting universe at late times with 
a < 0, d — 1. The physics of the early universe in this model is quite different from standard 
cosmology and requires a more detailed analysis including mattei^" however, it is a solid 



prediction that a{t) — )■ t for late times, i.e. the solution approaches a Milne universe. This is 
in remarkable good agreement with the basic observations including the type la supernovae 



data 



which are usually interpreted in terms of an accelerating universe. 



Milne universe. The Milne universe is nothing but (a quarter of) flat Minkowski space 



with fiat metric ds"^ 



-dr"^ + dr"^ + r'^dVt^ written in terms of the variables 
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r = tcosh(x), r = tsinh(x). 

Then this metric takes the form of a FRW metric with a{t) = t and k = —1, 

ds^ = -dt^ + a{t)\dx^ + smh\x)d^^)- (179) 

This is the metric in the above solution for late times, which is expected to be valid as long 
as the vacuum energy dominates the energy density due to matter. The solution is now 



easy to understand: it is an almost-trivial embedding 



C 



)10 



with a slight rotation in the 



early universe, like a twisted rubber sheet. Clearly such a solution is stabilized by large Ai. 

To obtain a reasonable theory of gravity, it remains to be shown that matter such as 
stars and galaxies lead to appropriate local perturbations of this cosmological solution, and 
reproduce at least Newtonian gravity. This indeed happens in this harmonic branch through 
an interesting mechanism, as we show next. 
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In particular the consistency with the CMB data cannot be rehably addressed at this point, see however 
for the simpUfied case of an exact Milne universe. The above refined solutions yield a big bounce and 

801, but the incorporation of matter is still missing. 



an early phase with power-law acceleration 
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7.2 Perturbations and Newtonian gravity from brane tension 

Now consider perturbations of the flat embedding M'^ C M"^ through the scalar fields 0*, 
leading to 

Qfiu = Vf^u + <9^0*<9i.0^ Sij = r]^y + h^y, (180) 

keeping only terms linear in h^^. Clearly we must keep contributions to 2nd order in 0* in 
order to get any non-trivial metric, hence Newtonian gravity will arise through a non-linear 
mechanism. We assume that G^y = g^^ as discussed above, and focus on static metrics 
g^y, corresponding to static and somewhat localized matter distributions. Thus consider the 
following localized excitation of the embedding 

t) = g{x)e^-' = g{x) ( ) , . = 1,2 (181) 

with very small u. This leads to the metric 

ds'^ = -(1 - g^) dt^ + (% + digdjg)dx'dx\ (182) 



which is static, because (181) is a standing wave rather than a traveling wave. 



7.2.1 Harmonic gravity bags 

Consider first vacuum excitations (l)Q{x,t) for p = 0. Neglecting corrections due to the 
induced gravity action, the equation of motion is simply 

□0^(x) = (183) 

which for spherical waves reduces to dr{r'^g')+u)'^r'^g{r) = 0. The unique localized spherically 
symmetric solution which is regular at the origin is 

, , sin(a;r) / cos(a;t) \ 

9o{r) = go^^, Mx)=9o{r){ ^,J, J ] (184) 

with radial wavelength given by 



sin(a;t) 



2ti 

UJ 



;i85) 



The effective metric (182) is 

ds^ = -(1 - u;2 ^(^^2) ^^2 ^ ^ {g'f)dr^ + r^rffi^ (186) 

which allows to read off the effective gravitational potential Uq seen by a static test particle: 
^oo = -(l + 2t/o), 2f/o(r) = ^uj'g{rf = -u'gl{^^^)\ (187) 



which satisfies 



T — y oo 

Uo{r) ~ <; r ; ^ ^ (188) 




Thus U{r) describes an attractive "gravity bag' with size L^^, decreasing as ^ for r > L^^. 
Due to the attractive gravitational force, matter will tend to accumulate inside these gravity 
bags. In particular, large clusters of matter such as galaxies will be embedded in such 
gravity bags. The essential point is that the matter within such a gravity bag will experience 
Newtonian gravity, due to a local deformation of the gravity bag. 



7.2.2 Perturbed gravity bags, Newtonian gravity and Poisson equations 

To understand the mechanism, consider first a spherically symmetric static mass density p 
around the origin within the radius vm- For r > tm, (f)'^{x) is again a solution of = 
and therefore must have the form 

0^ = ^7 r e-*, g(r)=go^ ^ ~ cos ^ + 189 

oor oor 

assuming ur <ti I. The phase shift 5 <^ 1 is due to the presence of matter at the origin, and 
is the key for obtaining Newtonian gravity. The effective metric ( 186 ) becomes 

goo = -1 + g^u' cos{26) + ^-^^^^^^ + 0((^)2) + 0(r) (190) 

and similarly for grr', the correction terms will be discussed later. This corresponds to a 
gravitational potential which for intermediate distances 

sin5<c<;r<l (191) 

is well approximated by a ^ potential with a constant shift, 

N 1 9 9 qnusm(26) 1 9 9 , . 

U{r) = --a;V ^ Uo- ^° ^ f/o ~ "2^^- (192) 

The phase shift 6 indeed turns out to be proportional to M, which gives Newtonian gravity. 
To understand this more generally, consider the e.o.m. (161) for 0* coupled to matter, 

□ 0' = ^ f^'dM', f"" = T^" - ^g^"" (193) 

replacing Dg ^ and using V^T^'^ = 0. Assuming p(x) = Too > and Tij pa in the 
presence of (non-relativistic) matteij^ this becomes 

A.{U(t)' = Mpdlcj)'. (194) 

The solution of this equation depends on the amplitude of the background "gravity bag" 
(poix) outside of the matter distribution. Since we are interested in the gravitational field 



^^This is not evident, as the simple embeddings below will generally not lead to Q^^ — 0. It might be 
justified either for more sophisticated embeddings, or - more interestingly - if A4 <C Apianck which is very 
appealing as we will see. 
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g(x) 




Figure 1: sketch of embedding function g{x) with short-scale perturbations and long-distance 
oscillations 



due to e.g. a star within a galaxy or some large-scale cosmic mass distribution, we make the 
ansatz 

(j)\x, t) = g{x)e'^\ g{x) = go{x) + 6g{x) (195) 

where \6g\ ^ go{x) ~ goiO) is varying on short scales according to p(x), while go{x) is slowly 
varying at the scale L^. Thus go{x) reflects the average mass distribution in the galaxy resp. 
a large cosmic structure, while Sg{x) is a small local perturbation due to e.g. a single star, 
as illustrated in figure [ij The separation g{x) = go{x) + 5g{x) corresponds to the splitting 
in (189). Now it makes sense to linearize in Sg{x); this is the crucial step. Then (194) gives 

{A + u')g = -8'Ku'^^g{x). (196) 

which for g{x) = go{x) + Sg{x) with 6g ^ go implies 

(A + uj^)6g = -87ru;2^((?o + Sg) ^ -8vru;'^ go- (197) 

The gravitational potential is then given by 

U{x) = Uo-uj^goSg{x) + 0{6g^) (198) 



(cf. (192)), which under the above assumptions satisfies 



(A + u')Uix) ^ -u'goiA + co')6g{x) = SttcoV,^. (199) 



This is the desired Poisson equation 



where we identify 



At/ ^ 47rG(p(x) + -^) = 47rGp(x) - Aeff (200) 
87r 



4 2 

G = 2^?o^ = -4^o^, (201) 
Aeff = -^GAt = 2Uoco^ (202) 
Uo = -\gW. (203) 
44 



Note that G is naturally small since \Uo\ < 1 and L^^ is large. The correction due to Af can 
be identified as effective vacuum energy resp. apparent negative cosmological constant Agfj. 
If p ^ this reduces to the usual Poisson equation of Newtonian gravity, AU = AnG p. 
Corrections to the Newtonian approximations will be discussed below. 

We conclude that localized matter p inside a gravity bag is subject to Newtonian gravity, 



with a dynamically determined gravitational "constant" G given by (201 ). For example, two 
stars or planets would lead to local perturbations 

g{x)= gQ + 5gi{x) + 5g2{x) (204) 

where 5gi are perturbations due to object pi. They both see the same gQ and cu, thus the 
same gravitational constant G, and Newtonian gravity is recovered on scales shorter than 
L^. However at very long scales and near the border of the galaxy resp. the gravity bag, 
the effective gravitational constant will vary. The role of the "extrinsic curvature" of the 



gravity bag and its brane tension will be clarified in section 7.4 using a different approach. 



Gravitational (non—) constant G. Since Li_j,go and therefore G are dynamical here, 
the question arises why G would be at least approximately the same in different parts of 
the universe. A possible explanation might be as follows: Galaxies are typically parts of 
large structures such as (super)-clusters and filaments, which should provide the dominant 
contribution to the brane embedding and hence to G. Given the homogeneity of the CMB 
background, it is plausible that similar scales for go{x) arise on these dominant cosmic 
structures. As a consistency check, one can show that u and go of gravity bags propagating 
in a fiat Milne resp. Minkowski background remain essentially constant in time p2]. 

Hence we obtain at least a crude approximation of the universe as we see it at late times, 
with far less fine-tuning than in the standard model. One may hope that a fully realistic 
picture will arise in a more sophisticated treatment. 



7.3 Beyond Newtonian gravity 

Now consider spherically symmetric mass distribution p(r) at the origin. A more careful 
analysis 22 shows that the time-component of the effective metric is given by 

^00 = - + 2f/o - 2^ + ^MGcV - ^Aeffr^ + (^)'^) + 0{r'). (205) 

We assume that the Newtonian potential due to M is smaller than the background potential 
Uq due to the harmonic bag, 

MG , , , , 

< \2Uo\. 206 

r 

Then 0{M^/r'^) term can be neglected, and the vacuum energy term (as well as the New- 
tonian potential) dominates the linear term, |Acfrr^| = 2|f/o|u;^r^ ^ ^^w^r^. Then goo 
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has approximately the form of a Schwarzschild-de Sitter metric with apparent negative 



cosmological constant Agg and a constant shift 84 



/ 2GM 1 . n\ , , 

(700 ^ - (l + 2f/o - - 3 Acfrr' j (207) 

assuming r < L^. The Newtonian term dominates the Aggr^ term provided the vacuum 
energy E^^{r) = ^^Af is smaller than M. We then obtain Newtonian gravity with potential 

U{r) ^ Uo-^. (208) 

The radial part of the effective metric turns out to be 

12GM 1 . GM 1 GM , o ^, ^ . 

g„^l + -Aesr^ -— 1 + cuV + O . 209 

3 r 9 r 2J\cE'r r 



7.4 



usmg a 



The factor | differs from general relativity, and will be confirmed in section 
different approach. This presents a challenge for the solar system constraints; however this 
might change in more sophisticated embeddings, and a more complete analysis is required 
before a reliable judgment can be given. 

The basic result is that Newtonian gravity arises at intermediate scales, with significant 
long-distance modifications. Notice that the precise form of the induced gravitational action 
was never used up to now, rather gravity arises through a deformation of the harmonic 
embedding which couples to T^^". Thus the mechanism is quite different from GR. 

Vacuum energy and cutoff. Inside the gravity bag, the vacuum energy > con- 



tributes a positive energy density to the gravitational potential (200) within the harmonic 
gravity bag, leading to an additional gravitational binding. However for very large distances 
r > L^, the harmonic behavior 



1 2 2^111^(1^'" 



Uir) - —^-^'ai^^ (210) 
dominates and leads to a screening of gravity, smoothly merging g^u{x) — > rj^y with the fiat 



metric of the Milne-like cosmology as discussed in section |7.1[ Recall that cosmology does 
not lead to the usual stringent constraints on the vacuum energy here. An illustrative plot 
of U{x) = —^u'^g{xY in comparison with the terms in (207) is given in figure [2j This clearly 



shows that the Newtonian potential dominates for small r, while the vacuum energy term 
takes over for larger r until the potential is cut off effectively at L^^. 

(Galactic) rotation curves. The non-relativistic orbital velocity v{r) around a central 
mass with the above metric is given for small distances r < L^^ by 



v(r) = VU^ = J2 — (1 + -— ) - -Aeffr2. (211) 
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Figure 3: Orbital velocity v{x) for a central point mass compared with Newtonian case: 
Newtonian domain, enhancement and cutoff, for u = 0.1, go = 1,6 = 0.1. 



This decreases like r~^^'^ as in Newtonian gravity as long as E^aci'^) < M, but for Eyac{f) ~ 
47rM it starts to increase linearly like v ~ a/IA^t until r L^. At that scale, the harmonic 
cutoff becomes effective, leading to a decreasing rotational velocity t> ~ ^ for large distances. 
An illustrative plot of v{r) compared with the Newtonian approximation is shown in figure 
|3j This idealized result hold only outside of the mass distribution, and will be modified 
e.g. in galaxies. One naturally obtains a slightly increasing rotation curve, which is indeed 
often observed. It remains to be seen whether this allows to explain the galactic rotation 
curves, providing an (partial?) alternative for dark matter. We only emphasize here that the 
qualitative behavior certainly goes in the right direction, and G = may differ somewhat 
from galaxy to galaxy. 



7.4 Local linearized effective action 

The above analysis requires a non-linear treatment of the large-scale gravity bag. We now 
provide a local, linearized description of the same mechanism Consider a local neighbor- 
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I am indebted to N. Arkani-Hamed for illuminating discussions and insights related to this section. 
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hood of a point p inside a gravity bag, using "free-falling" normal embedding coordinates 
9^J.4>o\p = 0, QiMulp = Vfj.u- Then the background gravity bag manifests itself in terms of a 
non- vanishing 2nd derivatives of the embedding' 
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d,d,<Pl = K;,, = + Oix'). (212) 

For example, a spherical harmonic gravity bag (j)o{x) = go cos{ut) ^^^^"^^ leads near its center 
to ^K^^x^^x" = —^gooj^it'^ + + y"^ + z^)). The advantage of this approach is that it 
allows to explore more general backgrounds, and indeed Ricci-flat metric perturbations are 
obtained for i^^^, ~ 6f^u. The linearized action valid near p is given by 

= J d^x {A{(P'D(f)^ 6ij - 87r0' f'^'d^d^cp^ 6ij) (213) 



using partial integration, recovering (193). We now add a small perturbation corresponding 
to a small mass distribution p, 

cf){x) = 0o(x) + 6(P{x), 000 = Snj^K^,. (214) 

which we assume to be static with T'^'^ = S^'^S^'^plx) for simplicity. We can then assume 
that 6(f){x) is static (more precisely d^Scp ^ -^^00)5 so that the effective action for S4>{x) can 
be written as 

^ j d^x {Al5(t)'U5(t)^ Sij - 16n5(j)' pK^^o ^ij) (215) 
omitting background terms. This is now reminiscent of scalar gravity where S(j) couples to 



p, cf. 85 . The equation of motion for 6(f)'' becomes 



AS(P' = Mf^'^K^,— ^ SnpKoo^, + (gravity bag), (216) 



cf. (197). The effective metric has the form 

g/Mu = + d^(f)ldy(f)l5ij + d^,(t)Qdy5(t)^ij + d^6(f)'dy(l)o6ij + df,5(f)'dj(l)^ij (217) 

The mixed goi components can be eliminated using a change of variables x'^' = x^ + = 
d^>(f)Q6(f>'6ij which gives 

g'^, = 7]^,-2KlJ<j)^S^j + d^S<j)'dM^5,,. (218) 



'^^KJ^^ is essentially the 2nd fundamental form. 
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Thus we recover the Newtonian potential U{x) = KQQ6(f>'6ij as in (198), which according to 
(216) satisfies 



1 



^^,y-^JJ- ^4 -47rGp + (gravity bag) 



(219) 



The gravitational constant G is determined by 



G 



'■oo-''-oo"u 



A! 



(220) 



in agreement with (201). In the case of a point mass M and a harmonic gravity bag, this 
gives 



5^00 



'1 



CM. 



g,, = 5ij{l + - ) + 0{5g'). 

O T 



{221] 



The factor \ agrees with our previous result (209), and stems from the harmonic form of the 



gravity bag. It is not hard to see ^81j that the leading Ricci-fiat metric perturbation 

1 



g' = 71,, - 2U{x)6, 



U{x) = ^K6(f) 



(222) 



of linearized GR is obtained for K,,^ = ^6,,K. For static 6(f), this gives for the Einstein 

tensor Qqq ~ SnGp and Q^P = 0, consistent with GR. In other words, a gravity bag with the 
form 

00 = cos(a;t)^^^^ ^ (1 - ^u\t^ + x^ + + z^)), 

with fc^ = 3uj^ would lead to a Ricci-fiat local deformation /i^j^; however it is not clear why 
such a background should arise. 

We conclude by pointing out that the effective gravitational constant G in the harmonic 
branch is governed by the extrinsic geometry A'^^ of a "gravity bag" , whose intrinsic geometry 
may be (almost) flat. This is an interesting mechanism which deserves to be studied in detail, 
and it remains to be seen whether a more sophisticated version will modify the detailed 
properties and perhaps reconcile it with general relativity. 



8 Conclusion 

The foundations of the matrix-model approach to (emergent) gravity are presented. The 
effective geometry of noncommutative space-time branes Aie C is identifled in the semi- 
classical limit, and the effective action for flelds, matter and geometry is given. Fluctuations 
of the matrices around such a background lead to flelds and matter on Aie, and the com- 
mutators in the matrix model become derivative operators which act on these flelds. Hence 
the noncommutative nature of Aig is essential for this mechanism, nevertheless 6^'^ is largely 
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hidden from low-energy physics and enters only through the effective metric. The matrix 
model action can be written in a covariant i.e. geometric way, leading to reasonable effective 
actions which can be studied further. In particular the D = 10 models, notably the IKKT 
model, allow to describe generic 4-dimensional geometries, and can be expected to define 
a well-behaved quantum theory. Moreover, the Einstein-Hilbert action are obtained either 
upon quantization or as an additional higher-order term in the matrix model. Therefore the 
matrix model should provide a realization of Einstein gravity or some closely related gravity 
theory. Emergent gravity from matrix models therefore becomes a serious candidate for a 
realistic theory of gravity at the quantum level. 

There are several aspects which make this framework very attractive from a theoretical 
point of view. First, its definition requires no classical-geometrical notions of geometry 
whatsoever. The geometry arises dynamically, which is very appealing from the point of 
view of quantum gravity and cosmology. Another fascinating aspect is that the matrix model 
framework leads naturally to a unified picture of gravitons and nonabelian gauge fields, which 
arise as abelian resp. nonabelian fluctuations of the basic matrices (covariant coordinates) 
around a geometrical background. The quantization around a such a background should be 
technically rather straightforward, similar to nonabelian gauge field theory. Nevertheless, 
more analytical and numerical work on the emergence and stability of 4-dimensional NC 



branes in matrix models is needed, particularly for the supersymmetric case 42 , 86 . 

Leaving aside its theoretical appeal, we can briefly summarize the physical aspects of 
the model as follows. 2 types of solutions have been identified, dubbed "harmonic branch" 
and "Einstein branch" . The harmonic branch essentially describes branes with tension. The 
most interesting feature is that it naturally predicts a (nearly-) fiat universe, resulting in 
luminosity curves e.g. for type la supernovae which are close to the observed ones (usually 
interpreted in terms of cosmic acceleration) without any fine-tuning. These and other long- 
distance modifications of gravity might also offer a mechanism for (partially?) explaining the 
galactic rotation curves, without requiring large amounts of dark matter. Moreover, fiat space 
remains to be a solution even in the preesnce of large vacuum energy, which is very interesting 
in the context of the cosmological constant problem. However, the solar system precision 
tests may be a problem in the harmonic branch. These should be obviously satisfied in the 
Einstein branch, which has not yet been studied in this framework in great detail; see [8l] 
for a realization of the Schwarzschild geometry. Generally speaking, the physical properties 
of emergent gravity are not yet sufficiently well understood. This applies in particular for 
the would-be U{1) gauge field which appear to behave as additional gravitational waves. 

We conclude by stressing that the gravity "theory" under consideration is not based on 
some theoretical expectations, but is simply the result of a careful analysis of the semi- 
classical limit of this type of matrix models. Taking simplicity as a guideline towards a more 
fundamental theory, the matrix model appears to be an extremely appealing candidate for 
a quantum theory of fundamental interactions including gravity, and certainly deserves a 
thorough investigation. 
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Appendix A: Some identities 

The following is a useful identity for Poisson tensors: 

= -d^^e^"'" -29''''p-'d,,p (223) 
where p = 6*"" is the symplectic volume, noting that 2p'^duP = duO^^' 9~^,. This imphes 

9^(p0n=O. (224) 
For our restricted class of metrics, this identity together with |G'^i/|^'^^ = pe'^ implies 

-e-"d,{pe''-'e^^'g^,Ax)) 



1 

— ( 
p 

e"'e'"''d,ie^^'''gMx)) (225) 



Appendix B: Identity for 0^^ 

To obtain a covariant equation for 9^'^, recall the following identity derived in ( [Ts], Appendix 
B) 

G-y^ix) v,9-^ = G^^ d,e-^ - T^^^e-;, - T^e-l 

which in (any) matrix coordinates can be written as 

G^\x) V^e-^ = Gp, 9"^ [e-^'^d^v + e"'^ d^x^ G^^d^d^x^ - e-^'G^'^g^^ d^a) - 

= Gp, 9"^ [e-^'^d^v + e"^ d^x'' (□gx'^ + T^d^x^) - e-^'G^^g^rj d^a) - r>'9;l 
= Gp, 9P^ [e-^'^d^T] + e-" d^x'' Ugx'^) 

+ (e-'^ Gp, 9>"^g,^)T^ - T^^9-J - G^^'(e-^ G,, 9P^g^,) d,a 
= Gp, 9^^" {e'^'^d^T] + e-" d^x" Dgx") + G^'< 9^^ d^a (226) 

which gives 

G^\x) V^(e'^0^;) = Gp, 9'^ {e-'^d^T] + dpx'' Dgx") (227) 
This is interesting, because it shows that 9'^'^ captures some of the extrinsic geometry of 



Ai C M'^ through the last term. In particular, the conservation law (107) implies 

d^x" Dgx" = (228) 
pi], which holds identically for g^i^ = Gp^. 
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